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Abstract 



We present exact calculations of flow polynomials F{G, q) for lattice strips of var- 
ious fixed widths Ly and arbitrarily great lengths L^, with several different boundary 
conditions. Square, honeycomb, and triangular lattice strips are considered. We in- 
troduce the notion of flows per face // in the infinite-length limit. We study the zeros 
of F{G, q) in the complex q plane and determine exactly the asymptotic accumulation 
sets of these zeros B in the infinite-length limit for the various families of strips. The 
function // is nonanalytic on this locus. The loci are found to be noncompact for many 
strip graphs with periodic (or twisted periodic) longitudinal boundary conditions, and 
compact for strips with free longitudinal boundary conditions. We also find the inter- 
esting feature that, aside from the trivial case Ly = 1, the maximal point, qcf, where 
B crosses the real axis, is universal on cyclic and Mobius strips of the square lattice for 
all widths for which we have calculated it and is equal to the asymptotic value qcf = 3 
for the infinite square lattice. 
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1 Introduction 



Consider a connected graph G — (V, E) with vertex set V and edge (bond) set E and an 
abehan group H of order o{H), represented as an additive group. For some general results, 
it will be necessary to allow the possibility of multiple edges joining a given pair of vertices, 
and loops (edges joining vertices to themselves), but these will usually not be present on the 
strip graphs of interest here. For definiteness, we let H be the additive group of integers mod 
q, 7jq. Denote the number of vertices and edges as n = \V\ and \E\. Assign an orientation to 
each of the edges in G and consider a mapping that assigns to each of these oriented edges 
a nonzero element in Z^. A flow on G is then defined as an assignment of this type that 
satisfies the condition that the flow into each vertex is equal to the flow outward from this 
vertex. Here the addition of two flows is deflned according to the abelian group H, so that, 
for H = Zq, the flow into each vertex is equal to the flow out of this vertex mod Z^. In a 
fluid or electric circuit analogy, this means that the (discretized) flow or electric current is 
conserved at each vertex mod q; there are no sinks or sources. Since the assignment of the 
zero element of to a given oriented edge is equivalent to the absence of the edge as far 
as the flow is concerned, a standard restriction is that an admissible flow must avoid zero 
flow numbers on any edge; this is termed a nowhere- zero g-flow. Henceforth, since all of our 
discussion will concern nowhere- zero g-flows, we shall take "g-flow" to mean "nowhere-zero 
g-flow". Some works discussing flow polynomials include [l]-[24]. 

An important problem in mathematical graph theory concerns the enumeration of g-flows 
on a given (connected) graph G. Tutte showed that there exists a polynomial in g, which we 
denote F{G,q), that equals this number of g-flows [3, 5]. The existence of this polynomial 
relies upon the fact that, for a given abelian group H, the number of flows on G depends 
only on the order of H, o{H) = g, not on any other structural properties of H. Clearly if 
G contains a bridge (isthmus), then there are no g-flows since the fluid or electric current 
flowing across the bridge has no way of returning. Thus, for any such edge that is a bridge, 
denoted e^, one has F{eij, g) = 0. The minimum (integer) value of g such that a bridgeless 
graph G admits a g-flow is called the flow number of G, denoted (t){G). 

In this paper we shall present exact calculations of flow polynomials for a number of lattice 
strips with various widths Ly and arbitrarily great lengths having several different types 
of boundary conditions. The lattices considered include the square (sg), honeycomb {he), and 
triangular (tri) lattices. The longitudinal (transverse) direction is taken as the horizontal, 
X (vertical, y) direction. We envision the strip of the triangular lattice as being formed by 
starting with a strip of the square lattice and then adding edges joining, say, the lower left 
and upper right vertices of each square. The strips of the honeycomb lattice are envisioned as 
brick lattices. We use the symbols FBCj, and PBC^^ for free and periodic transverse boundary 
conditions and FBCj,, PBCj,, and TPBCj; for free, periodic, and twisted periodic longitudinal 
boundary conditions. The term "twisted" means that the longitudinal ends of the strip are 
identifled with reversed orientation. These strip graphs can be embedded on surfaces with 
the following topologies: (i) (FBCy,FBCa;): free or open strip; (ii) (PBCy,FBCa;): cylindrical; 
(in) (FBCj/,PBC^): cyclic; (iv) (FBCy,TPBC^): Mobius; (v) (PBCy,PBC^): torus; and (vi) 
(PBCj„TPBC^): Klein bottle. 
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We shall introduce the notion of flows per face, // in the limit \V\ oo. From the 
basic definition one can generalize q from Z+ to R or, indeed, C. This generalization is 
necessary when one calculates the zeros of F{G, q) in the complex q plane. Using our exact 
calculations of F{G, q) for a variety of families of lattice strip graphs, we shall determine 
exactly the continuous accumulation set of these zeros of F{G,q) as \V\ ^ oo for each 
family of graphs G. This formal limit of the family of graphs G as \V\ ^ oo is denoted 
{G}. We shall call this accumulation set B or, when necessary to distinguish it from other 
accumulation sets, Bfi. This locus is of interest because the function fl is nonanalytic across 
B. We also study the approach of the fl functions for various infinite-length strips of the 
square lattice, as functions of q, to numerical values for the infinite square lattice. 

There are several motivations for this work. The fiow polynomial is of fundamental 
importance in graph theory. As we shall discuss below, it is a special case of the Tutte 
polynomial and encodes information on the connectivity of a graph. It is also of interest 
from the viewpoint of statistical mechanics because it is a special case of the g-state Potts 
model partition function. To our knowledge, there has not been any study of the function // 
that we introduce, and the associated regions of analyticity of // separated by the locus B 
for strips of regular graphs. In particular, the points where B crosses or intersects the real q 
axis are somewhat analogous to phase transitions in statistical mechanics, in the sense that 
the // function has different analytic forms on different sides of B. Associated with this, 
one finds qualitatively different flow behavior, as described by //, in the different regions 
separated by the locus B. While the ffow polynomial of a planar graph is equivalent to 
the chromatic polynomial of its dual graph, as we shall discuss further below, we find that 
the locus B = Bfi for the infinite-length hmit of a strip graph of a regular lattice with 
periodic (or twisted periodic) longitudinal boundary conditions is rather different from the 
corresponding locus of zeros Bw of the chromatic polynomial for the infinite-length limit of 
this strip. These differences are especially intriguing in view of a theorem that we shall give 
below, showing that in the limit of infinite width, i.e. for the two-dimensional lattice A, 
Bfi{A) — Bwi^*), where A* is the planar lattice dual to A. 

The fiow polynomial also appears in applied mathematics, engineering, business, and 
economics. Important questions in engineering and business concern the optimization of the 
fiow capacity through parts of networks and the study of where bottlenecks in fiows can 
occur because of limited connectivity. Although we shall not have occasion to make use of 
it, we mention the "max- fiow, min-cut" theorem of Ford and Fulkerson [6]. In the present 
work, we shall concentrate on intrinsic properties of fiow polynomials and their connections 
with statistical mechanics. 

We comment on the relation with previous calculations of chromatic polynomials (e.g., 
[25]-[62]). Given the connection (3.1.2), calculations of chromatic polynomials for families of 
planar lattice strip graphs immediately yield fiow polynomials for the dual graphs. We shall 
not recapitulate these calculations here since they are available in the literature. Chromatic 
polynomials for nonplanar graphs cannot be immediately related to fiow polynomials. 
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2 Connection with Tutte Polynomial and Potts Model 



In this section we review the connection between the flow polynomial and the Tutte poly- 
nomial or equivalently the Potts model partition function. Consider, as before, a connected 
graph G = {V^E). If an edge of G is a loop, e = e^, then the fluid (or electric current) conser- 
vation condition is automatically satisfied for any value of so that F(e£, q) = g — 1, where 
for (nowhere-zero) flows the zero element of H is excluded and there are thus q — 1 choices for 
the flow through the loop. Wc mention the standard notation that for an edge e E E, G/e 
is the graph obtained from G by contraction on e, i.e., deleting the edge e and identifying 
the vertices that it joins, and G — e is the graph obtained from G by deleting e. Then for an 
edge e that is not a bridge or a loop, the flow polynomial satisfies the deletion-contraction 
relation 

F{G, q) = F(G/e, q) - F{G - e,q) . (2.1) 

If G is a graph and is a field, then a general function g : G ^ K is a, Tutte- 

Grothendieck invariant if (1) if e G is a bridge, then g{G) = gbg{G/e), (2) if e e is a 
loop, then g{G) = gigiG/e), (3) if e e £^ is neither a bridge nor a loop, then 

g{G)^ag{G-e)+hg{G/e) a,b^O. (2.2) 

Prom the discussion above, it follows that F{G, q) is a Tutte-Grothendieck invariant. A 
useful property is that a Tutte- Grothcndieck invariant can be expressed in terms of the 
Tutte polynomial. A spanning subgraph G' = {V, E') of a graph G = {V, E) is a subgraph 
with the same vertex set and a subset E' Q E oi the edge set of G. The Tutte polynomial 
of a graph G is defined as 

T{G,x,y) = 5^ (x - l)HG')-kiG)^y _ i^c(G') (2.3) 

G'CG 

where k{G') and c(G') denote the number of components and linearly independent circuits 
of G', where, for an arbitrary graph G = {V, E), 

c{G) ^ \E\ - \V\ + k{G) (2.4) 

The quantity c(G') is also called the co-rank or nuUity of G', and the co-rank of the full graph 
G is called its cyclomatic number. The rank of a graph G is defined as r(G) = | — k{G), so 
that the first exponent can be written equivalently as k{G') — k{G) = r{G) —r{G'). Since we 
deal only with connected graphs here, k{G) = 1. (We follow the standard notational usage 
of X and y for the arguments of the Tutte polynomial and caution the reader not to confuse 
these with the x and y directions along the strip graphs.) 
A Tutte-Grothendieck invariant g is then given as 

g{G) = al^l-l^l+i6l^l-^T(G, ^. ^) . (2.5) 

b a 

In particular, for the fiow polynomial, in terms of the notation above, we have g^ — 0, 

gi = q — 1, a = —1, 6=1, whence 

F(G, q) = (-l)l^l-l^l+iT(G, X = 0, y = 1 - g) . (2.6) 



3 



Thus, the flow polynomial is a special case of the Tutte polynomial. Hence, we can use 
our previous exact calculations of Tutte polynomials and Potts model partition functions for 
various families of graphs [63]- [70] to derive flow polynomials for these graphs. 

From (2.3) and (2.6), it follows that the degree of F{G,q) as a polynomial in q is the 
cyclomatic number of G, 

degF(G,g)=c(G) (2.7) 

As is clear from the fact that the cyclomatic number c{G) gives the number of linearly 
independent circuits, it is closely related to the number of faces of G, f{G). In particular, 
for a planar graph G, using the Euler relation \ V\ — \E\ +f{G) = 2, we have f{G) = c{G) + l. 

Next, we recall the equivalence of the Tutte polynomial to the Potts model partition 
function. On a lattice, or more generally, a graph G, at temperature T, this model is deflned 
by the partition function [10] 

Z(G,g,^) = J]e-^^ (2.8) 

{cr«} 

with the Hamiltonian 

n^-jJ25^.a, (2.9) 

(ij) 

where CTj = 1, ...,q are the spin variables on each vertex i & G, and (ij) denotes pairs of 
adjacent vertices. We use the notation 

(where there should not be any confusion between the variable v and the edge set V or 
between the variable K and the number of components of a graph k{G)) so that the physical 
ranges are (i) > corresponding to oo > T > for the Potts ferromagnet, and (ii) 
— 1 < ti < 0, corresponding to < T < oo for the Potts antiferromagnet. 

As before, let G' — (V, E') be a spanning subgraph of G. Then Z(G, g, v) can be written 
as [71] 

Z{G,q,v)= q'^'^'^v"^''^ . (2.11) 

G'CG 

This formula enables one to generahze q from Z+ to IR or, indeed, C. Prom it one also 
directly infers the equivalence 

Z{G, q, v) = {x- l)^(^) {y - 1)"T(G, x, y) (2.12) 

where 

x=l + - (2.13) 

V 

and 

y^v + l (2.14) 

so that 

g = (x-l)(y-l) . (2.15) 
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Combining (2.6) and (2.12), one has the relation 

F{G,q)^{-l)\%--Z{G,q, -q) 



(2.16) 



Thus, the flow polynomial is a special case of the g-state Potts model given, up to the above 
pref actor, by the evaluation 

v^-q , i.e., K = ln(l - q) . (2.17) 

For the usual case of flows with q > 2, the condition v — —q corresponds to a complex- 
temperature regime for the Potts model. In eq. (2.17), one would thus have K — ln(g — 1) + 
(2£+ l)i7r, £ e Z. 

3 Some General Properties of Flow Polynomials 

3.1 Duality Relation with Chromatic Polynomial 

Before presenting our new results, we recall some basic information about flow polynomials 
that will be relevant for our work. We first discuss a useful connection with chromatic 
polynomials. The chromatic polynomial P{G, q) of a graph G counts the number of proper 
Q'-coloring of G, where a proper gr-coloring is a coloring of the vertices of the graph using q 
colors, subject to the condition that the colors assigned to adjacent vertices are different. 
The minimum (positive integral) value of q that allows one to carry out a proper g-coloring 
of G is called the chromatic number of G, x{G). The chromatic polynomial is given in terms 
of the Potts model partition function and Tutte polynomial as 

P(G, q) = Z{G, g, (-l)"+ig^(^)T(G, x = 1 - g, y = 0) . (3.1.1) 

Physically, the chromatic polynomial is the special case of the partition function for the 
Potts antiferromagnet at zero temperature, v — —1. Now let G be a planar graph and G* 
its planar dual graph. Then 

F{G,q) = q-'P{G\q) . (3.1.2) 

This follows directly from the expressions for the flow and chromatic polynomials in terms 
of the Tutte polynomial, together with the symmetry relation of a Tutte polynomial 

T{G,x,y)^T{G*,y,x) (3.1.3) 

for a planar graph G. The equality (3.1.2) will be important for our later discussion. Several 
corollaries follow from this equality. First, let G be a bridgeless planar graph and G* its 
planar dual, which thus has no loops. Then 

0(G) =X(G*) . (3.1.4) 

(The restriction that G has no bridges is made so that flows exist on G, and this is dual to 
the restriction that G* has no loops, so that proper g-colorings of G* exist.) 
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3.2 Some Results on Existence of g-Flows 



We mention here a few mathematical results on existence of g-flows that are relevant to 
our work. Let us denote the degree or coordination number of a vertex in a graph G (i.e., 
the number of edges connected to this vertex) as A. An elementary theorem states that a 
bridgeless graph admits a 2-fiow if and only if all of its vertex degrees are even. This can be 
understood as follows. After choosing an orientation for the edges, assign to each oriented 
edge the flow value 1. Since 1 = — 1 mod 2, and since 2n = mod 2, the fact that each vertex 
has even degree means that the flows into each vertex are mod 2. A related statement is 
that if a graph has flow number (p{G) = 2, then, since (for the nowhere-zero flows considered 
here) there is only one choice of flow number for each edge, namely 1, which is the same as 
-1 mod 2, it follows that F{G, g = 2) = 1. 

One of the most important theorems pertaining to flow polynomials for planar graphs 
is the 4-flow theorem, which states that every bridgeless planar graph has a 4-flow. This 
is equivalent to the celebrated 4-color theorem [72, 73], that every bridgeless planar graph 
has a coloring of faces with four colors such that any two faces that are adjacent across a 
given edge have different colors. An equivalent expression of the theorem for proper vertex 
colorings is that every loopless planar graph has a proper 4-coloring. One can inquire about 
the existence of flows for arbitrary, not necessarily planar, graphs. Jaeger proved that every 
bridgeless graph (not necessarily planar) has an 8-flow [7] , and subsequently Seymour proved 
the stronger result that every bridgeless graph has a 6-flow [8] . An outstanding conjecture, 
due to Tutte [3] , is that every bridgeless graph has a 5-flow. These general results provide a 
useful background for the specific flow numbers that we shall obtain for various planar and 
nonplanar families of graphs. 

3.3 Structural Properties for General Strip Graphs 

Since a flow polynomial is a special case of a Tutte polynomial or equivalently, a Potts model 
partition function, one can obtain several properties of flow polynomials from corresponding 
properties of the latter two polynomials. A general form for the Tutte polynomial for the 
strip graphs considered here, or more generally, for recursively defined families of graphs Gm 
comprised of m repeated subunits (e.g. the columns of squares of height Ly vertices that are 
repeated = m times to form strip of a regular lattice of width Ly and length (denoted 
Ly X Lx) with some specified boundary conditions), is [63] 

T{Gm,x,y) = -— r E cg,,(At,gj)™ (3-3.1) 

where the terms \T,G,j^ the coefficients Cgj, and the total number Nt^g.x depend on G 
through the type of lattice, its width, Ly, and the boundary conditions, but not on the 
length. Equivalently, 

Z{Gm,q,v)^ J2 CGA^z,G,jr (3.3.2) 
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where 

Nz,G,x = Nt,ga ■ (3.3.3) 

It follows from (2.6) and (3.3.1) that the flow polynomial for recursive families of graphs, 
comprised of m repeated subgraph units, has the general form 

FiG^,q)= J2 ^GA^F,G,r (3.3.4) 

where again the terms XF,G,j, the coefficients cqj, and the total number of terms A^f,g,a 
depend on G through the type of lattice, its width, Ly, and the boundary conditions, but 
not on the length. This is related to (3.3.1) as follows. The prefactor in (2.6), (— l)l-^l-l^l+i 
is of the form (— 1)^^™+^, where c is an even or odd integer. The single factor (—1) cancels 
the (—1) factor resulting from the x = evaluation of the overall prefactor l/{x — 1) in 
(3.3.1). For the cyclic lattice strips considered here, both \V\ and \E\ are integer multiples 
of Lj. = m, so their difference is of the stated form, cm. Thus, 

AF,G,,(g) = {-irXT,G,j{x = 0,y=l-q) (3.3.5) 

for the subset of the At.gj's that are nonzero when evaluated at a; = 0. 

For a given type of strip graph, the sum of the coefficients in (3.3.4) is denoted 

Cf,g = Yl ^G,3 ■ (3.3.6) 



3.4 The Function // and Associated Locus B 

Given the structural property (3.3.4), it is natural, in the limit | V| — oo, to deflne a function 
specifying the number of flows per face, 

//({G},g)= hm F(G,g)i/^(«) (3.4.1) 

\V\-^oo 

For almost all of the strip graphs considered here, |V^| — > oo <(=^ f{G) oo. An exception 
is the circuit graph, for which F{G) — 2, independent of \V\. 
For reference, the reduced free energy is 

/({G}, q,v)^ hm HZ{G, q, t;)^/l^l] (3.4.2) 
|y|— >oo 

a limiting function obtained from the Tutte polynomial may be written as 

r({G},x,y)= hm T(G,a;,|/)VI^I (3.4.3) 
I y |— »oo 

and the ground state degeneracy per vertex of the gr-state Potts antiferromagnet is 

W{{G},q)= lim P(G,g)^/l^l . (3.4.4) 

|y|— >oo 
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Since q) is a polynomial in q of degree c{G), it follows that with this definition, 



■in. = (3.4.5) 

q-^oo q 

This is the analogue to the property that 

lim l^iMlil = 1 . (3.4.6) 

q^oo q 

Following our earlier nomenclature [32], we denote a term A in (3.3.4) as leading (= 
dominant) for a given value of q if it has a magnitude greater than or equal to the magnitude 
of other A's evaluated at this value of q. In the limit \V\ — > oo the leading A in F{Gm: Q.) 
determines the function fl{{G}, q). The continuous locus B where fl{{G}, q) is nonanalytic 
thus occurs where there is a switching of dominant A's in F and is the solution of the 
equation of degeneracy in magnitude of these dominant A's. Since a zero in F{G, q) requires a 
cancellation between dominant A's, it is clear that as —>■ oo, the continuous accumulation 
set B of the zeros oi F{G,q) forms the boundary curves across which this switching occurs 
and hence across which fl{{G}) is nonanalytic. Depending on the family of strip graphs, 
the locus B may or may not separate the complex q plane into different regions and may 
or may not cross the real q axis. If it docs, we denote the maximal point where it crosses 
this axis as qcf{{G}). We also denote the region including the positive real q axis extending 
down from g = oo (and terminating at qcf{{G}) if the latter point exists) as region Ri. This 
region Ri is understood to include the maximal area to which one can analytically continue 
fl({G}) from the large-g positive real axis. For families of strip graphs where B separates 
the q plane into different regions, the function fl{{G}) has different analytic forms in the 
different regions. For a strip graph G, in region Ri, 

fli{G},q)^XR^,aom. (3.4.7) 

where Xm^dom. denotes the dominant A in region Ri. 

Just as was true for these functions /, r, and W, there are two subtleties in the definition 
(3.4.1): (i) which of the c{G) roots to take in the equation, and (ii) the fact that at certain 
values of q, one has the noncommutativity of limits 

lim lim F(G',g)^/-^(^) 7^ lim lim F{G,qy/f^^K (3.4.8) 

|y|^ooq— »9s q^qs \V\-*oo 

We have discussed these before in the context of the definitions of W and / [32, 63, 47]. 
Concerning item (i), for sufficiently large real q, in the region Ri, F{G,q) is real and posi- 
tive, so one chooses the canonical root for //({G}), which is real and positive. However, if 
B separates the q plane into different regions, then in the latter regions, there is no canon- 
ical choice of phase that one can make for the root (3.4.1) and hence only the magnitude 
\fl{{G}, q) \ can be determined unambiguously. For item (ii) the noncommutativity typically 
occurs at special values qs = 1, 4^{G) — 1. To avoid isolated discontinuities that would oth- 
erwise result where fl vanishes, we shall, following our earlier practice for W and / [32, 63] 
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and adopt, for the cases where such noncommutativity occurs, the second order of hmits in 
(3.4.8) 

fl{{G},qs)= Mm Urn ^(Cg)^/^^^) . (3.4.9) 

From the duahty relation (3.1.2) it follows that if G is a planar graph and G* is its 
planar dual, then, taking into account that the number of faces of G is equal to the number 
of vertices of G*, one has 

fl{{G},q)^W{{G*},q) . (3.4.10) 

Note that if we had defined fl{{G},q) as lim\v\^oo F'{G,qy^^^\ we would have obtained a 
different relation fl{{G}, q) = W{{G*}, q^, where p is a geometric factor depending on the 
particular graphs G. Since the infinite square lattice is self-dual, it follows that 

fl{sq,q) = W{sq,q) . (3.4.11) 

Since the infinite honeycomb and triangular lattices are the duals of each other, we also have 

fl{tri, q) = W{hc, q) , fl{hc, q) = W(tri, q) . (3.4.12) 

The dual of the kagome lattice is the diced lattice [74, 75], and hence 

fl{kag, q) = W{diced, q) , fl{diced, q) = W{kag, q) . (3.4.13) 

We next give some results for the locus B. Because of the fact that these A's arc degenerate 
in magnitude on B, it follows that fl{{G}, q) is nonanalytic but continuous across B. A basic 
property is that B{{G}) is invariant under complex conjugation, 

B{{G}) = B{{G}r . (3.4.14) 

This is a consequence of the property that the coefficients of each term in the flow polynomial 
are real and hence the zeros of the flow polynomial arc invariant under complex conjugation. 
Hence, the same property holds for their asymptotic accumulation set as |y| cx3. 

We next point out another consequence of (3.1.2) and (3.4.10) for planar graphs. For this 
purpose, we must append subscripts to distinguish Bfi and Bw, the respective continuous 
accumulation sets of zeros of the flow and chromatic polynomials in the limit \V\ — > oo. 
Then for the |y| — > oo limits of planar graphs {G}, 

BM{G}) = Bw{{G*}) 

Several corollaries follow. Since the inflnite square lattice is self-dual, 

Bfi{sq) = Bwisq) 

Since the inflnite triangular and honeycomb lattices are dual to each other, we have 

Bfi{tri) = Bw{hc) , Bfi{hc) = Bw{tri) (3.4.17) 



(3.4.15) 



(3.4.16) 
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Since the kagome and diced lattices ([74, 75]) are dual to each other, 



Bfi{kag) — Bw{diced) , Bfi{diced) — Bwikag) 



(3.4.18) 



Now consider the lattice strip graphs G[Ly x Lx,cyc] and G[Ly x Lx,Mb.]. We have 
shown earher [45, 47, 52] that the locus Bw is the same for the respective infinite-length 
hmits of both of these strip graphs. Wc also find that the same property holds for the Bfi 
computed here, but note that this does not follow from our previous results for Byy since the 
Mobius strips are nonplanar, and hence one cannot apply the duality relation (3.1.2) and 
the consequent equalities (3.4.10) and (3.4.15). We also find that for the families of strip 
graphs that we have calculated, Bfi is the same for boundary conditions corresponding to 
embedding on surfaces with torus and Klein bottle topologies. 

3.5 Noncompactness of B for Classes of Strip Graphs 

An important feature of the loci B = Bfi for the infinite-length limits of the families of 
lattice strips with periodic (or twisted periodic) longitudinal boundary conditions is that 
many of these are noncompact in the q plane, containing curves that extend infinitely far 
away from the origin. A fortiori, this means that for this families of graphs, as L^, — > cxd, 
there is no upper bound on the magnitudes |g| of the zeros of the flow polynomial. The 
noncompactness of Bfi for many strips of regular lattices is quite different from the loci 
Bw for the corresponding strips, which, as we showed in our previous works on chromatic 
polynomials, were compact in the q plane [32]- [61]. Sokal has proved that for an arbitrary 
graph G, any zero of the chromatic polynomial satisfies the following upper bound [48] 



where here A denotes the maximal vertex degree in G. (Indeed, our explicit calculations 
for a variety of families of strip graphs yielded loci Bw on which the values of max(|g|) 
were considerably less than the above-mentioned upper bounds for the respective strips.) 
As was discussed in [37, 38, 39], the locus Bw is noncompact if and only if the degeneracy 
condition in magnitude of two or more different dominant A's can be satisfied for arbitrarily 
large Clearly, the same condition holds for Bfi. One of us (RS), with Tsai, studied the 
conditions under which this could occur and constructed several general classes of families 
of graphs that were designed to satisfy this condition and thereby yield noncompact loci Bw 
[32, 37, 38, 39] (sec also [30, 48]). As noted in [37, 38], the condition that \iTa\v\^oo^ = oo 
is a necessary but not sufficient condition for Bw to be noncompact; an example of a graph 
with a maximal vertex degree that goes to infinity as |y | ^ oo but with a compact Bw is the 
wheel graph Ki + Cm- Here, Cm denotes the circuit graph, Kp denotes the complete graph 
on p vertices, defined as the graph each of whose vertices is adjacent to every other vertex, 
and the "join" G + H is defined as the graph formed by connecting each of the vertices of G 
to each of the vertices of H. (The B in this case is simply the unit circle |q' — 2| = 1.) 

Using the duality relation (3.1.2) for planar graphs, one can understand why Bfi is 
noncompact for the simplest nontrivial cyclic strip of the square lattice, namely for width 
Ly — 2. To do this, we observe that the planar dual of this 2 x cyclic strip is the circuit 




(3.5.1) 
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graph with vertices, augmented by having an extra vertex connected by edges to each 
of the vertices on the upper side of the strip and, separately, an extra vertex connected by 
edges to each of the vertices on the lower side of the strip. This is a special case of one of the 
families of graphs that were shown in [37, 38] to yield, in the oo hmit, noncompact loci 

i3, namely the family {Kp)^ + G,., where, h means that h of the edges in the complete graph 
Kp are cut, and G^. denotes an arbitrary r- vertex graph. We refer the reader to section 11 of 
[37] and section 2 of [38] for further details. Specifically, we find that the planar dual graph 
to the Ly — 2 cyclic stripf of the square lattice is the special case of {Kp)ij + Gr with p — 2, 
6 = 1, signifying the cutting of the single edge in the K2, and Gr = Cl^, the circuit graph. 
Using the duality relation (3.4.15), we therefore prove that Bfi for the ^ 00 limit of the 
cyclic Ly X Lx strip is noncompact, extending infinitely far from the origin in the q plane. 
Similar arguments give insight into the noncompactness of B for the wider strips considered 
here. Because of the noncompactness of Bfi in the q plane and the property that we find 
that this locus does not pass through q — 0, it will often be convenient to display it in the 
plane of the inverse variable 

u^^ (3.5.2) 

where it is compact. 

We note that not all strip graphs with periodic longitudinal boundary conditions lead to 

noncompact Bfi. An example of a family with compact Bfi is provided by the cyclic self-dual 
strips of the square lattice Gu{Ly x L,j.) that we studied in [57, 68]. Since these are self-dual 
planar graphs, one can immediately use (3.4.15) to infer that the loci Bfi are identical to the 
loci Bw given in [57, 68], which are compact. 

In contrast to the case with many families of strip graphs with periodic (or twisted 
periodic) longitudinal boundary conditions, we find that B is compact for the ^ 00 limit 
of strips with free longitudinal boundary conditions, as will be illustrated with specific exact 
solutions below. This feature is similar to our previous findings for Bw- Henceforth, where 
no confusion will result, we shall drop the subscript // on Bfi. 

4 General Structural Results for Cyclic Strips of the 
Square and Honeycomb Lattices 

In [51] it was shown that for cyclic and Mobius strips of the square lattice of fixed width Ly 
and arbitrary length L^ (and also for cyclic strips of the triangular lattice) the coefficients Cj 
in the Tutte polynomial are polynomials in q with the property that for each degree d there 
is a unique polynomial, denoted c^'^\ Further, this was shown to be 

c^'^ = u,M'"m = Y.{-^y {^"^ r (4.1) 

i=o \ J / 

where Un{x) is the Chebyshev polynomial of the second kind. A number of properties of 
these coefficients were derived in [51]. We list below the specific c'^'^^'s that will be needed 
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here: 

c(o) = 1 , c(i) = g - 1 , c(2) = g2 - 3g + 1 , (4.2) 

c(3) = g3 _ ^ _ ^ _ (^4 3^) 

Thus, the terms XT,Ly,j that occur in (3.3.1) can be classified into sets, with the \T,Ly,j{l-i v) 
in the (i'th set being defined as those terms with coefficient c^*^'. In Ref. [51] the numbers 
of such terms, denoted nT{Ly, d), were calculated. Labelling the eigenvalues with coefficient 
c^'^^ as \T,Ly,d,j with 1 < J < nT{Ly,d), the Tutte polynomial for a cyclic strip graph of 
length Lx — m can be written in the form [51] 

Ly nT(Ly,d) 

T{G,[LyXm;cyc.],x,y)^-—Y,c^'^ E (At,g.,l„.j)"^ ■ (4.4) 

d=0 j=l 

For the Mobius strip of the square lattice the coefficients may be either +c^^^ or — c*^*^) and 
the terms are accordingly labelled as \T,Ly,d,±,j, where 1 < j < nT{Ly,d,±). We have [51], 
using the notation TPBCy for twisted periodic b.c. in the longitudinal direction, 

nT{Ly,d,+) riT{Ly,d,-) 



■,3 

d=0 ' j=l j=l 



r(.g[L,xm;M6.],x,y) = ^£c('^)[ ^ (Aw,+,,)'" - ^ (^W,- 



(4.5) 



2 for even Ly 



where 

C^max = { ■ (4.6) 
[ ^ for odd Ly 

The number ririGs, Ly, d) of A's with a given coefficient c'^'^^ is [51] 

nT(G„ L„ c?) = ^^^^^^1^^ (^^^^ ^ for = sq, trt, he; < d < Ly (4.7) 

and zero otherwise. The total number NT,Ly,\ of different terms \T,Ly,j in eq. (3.3.1) for 
cychc (or Mobius) strips Gs of the square, triangular, and honeycomb lattices is [51] 

Ly 

NT,G,,r.y,x^J2''T{Gs,Ly,d) (4.8) 



d=0 



which was calculated to be [51, 66] 

NT,G.,Ly,>. = (^^^j) . (4.9) 

For arbitrary Ly, eq. (4.7) shows that there is a unique \T,Ly,d corresponding to the coefficient 
(Ad) q£ highest degree, d = Ly, and this term is 

XT,Gs,Ly,d=Ly = 1 ■ (4.10) 
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Hence, 

^F,sq,Ly,d=Ly = ^F,hc,Ly,d=Ly = ( — 1)^^'^^ (4-11) 

and 

^F,tri,Ly,d=Ly = 1 (4-12) 

(independent of Ly). Since this eigenvalue is unique, it is not necessary to append another 
index, as with the other A's, and we avoid this for simphcity. 

From (4.4), it follows that the flow polynomial for a cyclic strip of the Gg-type lattice 
also has the same type of form, i.e., 

Ly np(Gs,Ly,d) 

F{Gs[Ly X m; cyc.],x, y) = c^"^ Yl i>^F,GM^r ■ (4-13) 

d=0 j=l 

Let us denote the number of XF,Gs,Ly,d,j with coefficient c^'^^ as npiGg, Ly,d). Then for this 
strip 

Ly 

CF,G,,Ly = J2 ^^(^- ' ^)^^'^ ■ (4-14) 

d=0 

We now concentrate on cyclic strips of the square and honeycomb lattice. Aside from de- 
generate cases, the dual of the cychc strip of the square lattice of width Ly > 2 and length 
Lx is a strip of the same lattice with width Ly — 1 and length L^ with all of the vertices 
along the upper and lower sides connected to two respective external vertices. Now consider 
proper colorings of this dual graph. The total number of proper g-colorings of a transverse 
slice is q{q — It is elementary to show that this also holds for Ly = 1. From (3.1.2), it 
follows that 

CF,sq,Ly = (g - 1)^^ . (4.15) 
A similar argument shows that (4.15) also holds for cychc strips of the honeycomb lattice. 

We now derive the following structural theorem. 

Theorem 4.1 Consider the flow polynomial for a cyclic strip graph Gg of the square 
or honeycomb lattice of fixed width Ly and arbitrarily great length L^. For brevity, set 
npiGsi Ly, d) = np{Ly, d). The np^Ly, d), d = 0,1, ..Ly are determined as follows. One has 

npiLy, d) = for d > Ly , (4.16) 

nF{Ly,Ly)^l (4.17) 
71^(1, 0) = (4.18) 
with all other numbers npiLy, d) being determined by the two recursion relations 

nF{Ly + l,0)=nF{Ly,l) (4.19) 

and 

nF{Ly+l,d) — nF{Ly,d—l)+nF{Ly,d)+nF{Ly,d+l) for Ly > 1 and 1 < d < Ly + 1 . 

(4.20) 
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Table 1: Tabic of numbers np{Ly, d) and their sums, NF,Ly,\ for cyclic strips of the square and honeycomb 
lattices. Blank entries are zero. 



Ly i d ^ 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




1 





1 




















1 


2 


1 


1 


1 


















3 


3 


1 


3 


2 


1 
















7 


4 


3 


6 


6 


3 


1 














19 


5 


6 


15 


15 


10 


4 


1 












51 


6 


15 


36 


40 


29 


15 


5 


1 










141 


7 


36 


91 


105 


84 


49 


21 


6 


1 








393 


8 


91 


232 


280 


238 


154 


76 


28 


7 


1 






1107 


9 


232 


603 


750 


672 


468 


258 


111 


36 


8 


1 




3139 


10 


603 


1585 


2025 


1890 


1398 


837 


405 


155 


45 


9 


1 


8953 



Proof. We substitute for c^''^ from eq. (4.1) in eq. (4.15). We obtain another equation 

by differentiating this with respect to q once; another by differentiating twice, and so forth 
up to Ly-fold differentiations. This yields Ly + 1 linear equations in the Ly + 1 unknowns, 
np{Ly, d), d = 0,1, Ly. We solve this set of equations to get the np^Ly, d). □ 

Corollciry 4.1 

npiLy, Ly-l)^Ly-l. (4.21) 

We note that the recursion relations (4.19), (4.20) for the numbers nF{Ly,d) for cyclic 
strips of the square and honeycomb lattices are the same as the recursion relations that we 
derived earlier in eqs. (3.14) and (3.15) of [51] for the corresponding numbers np{Ly,d) 
for cyclic strips of the square and triangular lattice. Thus, the differences in the values of 
nF{Ly,d) for cyclic strips of the square and honeycomb lattice and the np{Ly,d) for cyclic 
strips of the square and triangular lattice are due to the different initial values of these 
quantities, i.e. (4.18) in the former case and np(l, 0) = np(l, 1) = 1 in the latter case. 

Corollairy 4.2 

[(i, + l)/2] . , i\ /r ■ i\ 



Proof. This follows immediately from the solution to our general recursion relations 
(4.19), (4.20). □. 
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Summing the npiLy, d) over d for a given strip with Ly, we obtain 




(4.23) 



where [x] is integer part of x. This total number also applies to Mobius strips of the square 
and honeycomb lattice. 

A generating function for the NF,Ly,x is [76, 77] 



From this, it follows that the number Np^i^^x grows exponentially fast with the width Ly of 
the cyclic strip of the square or honeycomb lattice, with the leading asymptotic behavior 



Analogous structural results can also be given for cychc strips of the triangular lattice. We 
proceed to present exact calculations of flow polynomials for a number of different lattice 
strips. 

5 Cyclic and Mobius Strips of the Square Lattice 

5.1 General 

We denote the cyclic family as sq{Ly, L^, FBCx, PBCy) = sq{Ly, L^, eye.) and the Mobius 
family as sq{Ly, L^, FBCx, TPBCy) = sq{Ly, L^, Mb.). The cyclic and Mobius strips of the 
square lattice of width Ly and length L^ have \V\ = L^Ly vertices and \E\ = Lx{2Ly — 1), so 
that the prefactor in (2.6) is (^—l)^='(^y~'^)+^ and hence the \F,sq,Ly,d,j are given by the subset 
of the \T,sq,Ly,d,i^ ilisA, are nonzero when evaluated for x = and y — 1 — q, multiplied by 
the prefactor (— 1)^''"^. Among these famihes of graphs, the lowest case, = 1, is just the 
circuit graph Cm, and the result is elementary; the only \F,sq,j is unity, and F{C„i, q) = 1 — ^ 
(independent of m). This polynomial has only a single zero, at g = 1. With the definition 
(3.4.1), we obtain fl{q) = \/q — 1. Henceforth, for brevity of notation, where no confusion 
will result, we shall omit the F in \F,sq,Ly,d,j and write this simply as Xsq,Ly,d,j, and similarly 
for other lattice types. 

5.2 Ly = 2 Cyclic and Mobius Strips of the Square Lattice 

Using the duality relation (3.1.2), one has 




,y — 



(4.24) 



A^F,L„A ~ L-^'^ 3^^ as ^ oo . 



(4.25) 



F{sq[2 xLx = m, eye], q) = q-^P{{K2)i + C^, q) 



(5.2.1) 
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where Kp is the complete graph, G + H is the join of the graphs G and {Kp)i, denotes the 
graph obtained by removing b edges from Kp, and Cm is the circuit graph with m vertices. 
In [32, 37, 38], one of us (RS), with Tsai, calculated the chromatic polynomial for the family 
(^2)1 + Cm- This immediately gives the flow polynomial, via eq. (5.2.1), with 

A.,,2,0,1 = ? - 2 (5.2.2) 

Kq,2,i,i = g - 3 (5.2.3) 

A.,,2,2 = -1 . (5.2.4) 

Hence, we find 

F{sq[2 X m, eye], q) = {q - 2)"" + c«(g - 3)" + c(')(-l)™ . (5.2.5) 

Note that nF{2,0) = nF{2,l) = ni?(2,2) = 1 and A^f,2,a — 3, in accord with our general 
structural formulas given above. 

By using results that were obtained in [63] for the Tutte polynomial, we obtain 

F{sq[2 X m, Mb.],q) = (g - 2)"^ + c^^\q - 3)"^ - (-I)'" . (5.2.6) 

Thus, 

nF{sq,2,0,+) = l, (sg, 2, 0, -) = 1, (sg, 2, 1, +) = 1, (sg, 2, 1, -) - (5.2.7) 

with npisq, 2, d,±) — for d>2. 

We observe that F{sq[Ly = 2, = m, eye], q) and F{sq[Ly = 2,Lx — m, Mh.], q) always 
have the factors (g — l)(g — 2). For m > 3 odd, F{sq[Ly = 2, L^, = m, eye.], q) also has the 
factor (g — 3), while for m even, F{sq[Ly = 2,Lx = m, Mb.], q) also has the factor (g — 3). 
These results show that 

{3 if m is even 
(5.2.8) 
4 if m > 3 is odd 

(for m—1, this flow polynomial vanishes since the graph contains a bridge). Further, 

{4 if m is even 
(5.2.9) 
3 if m is odd 

We recall that a graph G is /c-critical if P{G, q = x(G^)) = k\, where the chromatic number 
of G, x{G), was defined above. Somewhat analogously to the concept of /c-critical graphs for 
proper vertex coloring and chromatic polynomials, one may ask whether for a graph G with 
(p{G) = /c it is true that F{G, q = (p{G)) has a fixed value depending on q, i.e. for q = 4>{G), 
the number of flows on G is a fixed number rather than growing (exponentially) with \V\. 
Using our results we can answer the question for these cyclic and Mobius strip graphs. First, 
we recall the elementary observation that if 0(G) = 2, then F{G, 2) — 1. From eqs. (5.2.8) 
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and (5.2.9) we know that the flow number is 3 for the Ly = 2 cyclic and Mobius strips of 
the square lattice with even and odd L^. = m, respectively, and we find 

F{sq[2 X m, eye], q = 3) = 2 for even m >2 (5.2.10) 

and 

F{sq[2xm,Mb.],q = 3) = 2 for odd m > 1 . (5.2.11) 

In contrast, it is readily verified from our exact results that for the Ly = 2 cyclic and Mobius 
strips of the square lattice with odd and even — m, respectively, for which (j) — 4, the 
evaluation of the respective flow polynomials at q — A yields values that grow (exponentially) 
with the length of the strip. 

One way to prove eqs. (5.2.10) and (5.2.11) is to use our explicit calculations of the flow 
polynomials for these families of graphs. Another way is to write down the actual flows. For 
this purpose, consider a given square on the strip. For q — 3 there are two circular flows 
on this square, namely those with a flow number 1 assigned to each edge, going in (i) a 
clockwise or (ii) a counterclockwise manner. Clearly, a flow with flow number 1 assigned to 
each edge, going in a clockwise manner is equivalent to a flow with flow number 2 = — 1 mod 
3 assigned to each edge going in a counterclockwise manner. Consider the flow of type (i) 
on this square. In order to satisfy the flow conservation condition at the right-hand upper 
and lower vertices, it is necessary that the flow on the neighboring square to the right be of 
type (ii). Similarly, if the flow on the given square is of type (ii), then it is necessary that 
the flow on the neighboring square to the right be of type (i). Continuing in this manner 
along the cyclic strip, one finds a consistent set of choices if m is even, but not if m is odd. 
This proves the first part of the corollary. A similar proof with obvious changes works for 
the second part. 

We remark that the constructive proofs given above establish that ior q — 3 the total 
flows can be decomposed into superpositions of circulations around each square. Since there 
is a 1-1 correspondence with the flows around a face of a graph and the face- colorings of 
the graph [3], one sees that these corollaries are equivalent to an analogous condition for the 
face-coloring of the strip graph, or equivalently, the vertex coloring of the dual graph. 

In accordance with our general discussion given above, the locus B for the ^ oo limit 
of the cyclic and Mobius Ly = 2 strips of the square lattice is noncompact in the q plane. 
From eq. (5.2.1), it follows that this locus is identical to the one that was determined earlier 
for the m — > oo limit of the family P{{K2)i + Cm, q) in [32, 37, 38] and shown in the q plane 
shown in Fig. 2 of [32], or equivalently, to the locus B in the u — 1/q plane shown in Fig. 
1 of [38]. It divides the q plane into three regions, Ri, i = 1,2,3. Regions -Ri, i?2, and 
contain the respective real intervals 3 < q < oo; 2 < q < 3, and q <2. Region R2 is bounded 
on the left by the arc of the circle 

q^3 + e'\ for 
and on the right by the arc of the circle 

q = 2 + e''^ , for 



2Tvi ^ Am 

— <e < — 
3 - - 3 



3 - - 3 



(5.2.12) 



(5.2.13) 
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Thus, 

ge/ = 3 for {G} = sq, 2 x oo, cyc./Mb. (5.2.14) 
The two arcs (5.2.12) and (5.2.13) meet at the complex-conjugate points 

quq; = l±^ = 2 + e^-f' (5.2.15) 

which are triple points on the locus B. Curves on B extend upwards from qt and downwards 
from g* separating region Ri from R^. The curves on B pass through the origin u = 
verticaUy. This foUows because at u = there are two A's which are leading and are 
degenerate in magnitude. 

In region Ri, Xsq,2,o,i is dominant, so that 

fl{sq[2x oo, eye], q) ^ fl{sq[2x oo, Mb], q) ^ q -2 for g e i?i . (5.2.16) 

In region R2, Kq,2,2 is dominant, so that 

\fl{sq[2x 00, cyc./Mb.], q)\ = l for g e i?2 ■ (5.2.17) 

In region R^, Xsq,2,i,i is dominant, so that 

|//(sg[2 X oo,q/c./M6.],g)| = |g-3| for g e i?3 . (5.2.18) 



5.3 Ly = 3 Cyclic and Mobius Strips of the Square Lattice 



Our general structural results yield np{3,0) = 1, ^^^(S, 1) = 3, np{3,2) = 2, ni?(3,3) = 1 
with the total number Np^^^x = 7. For this Ly = 3 case we find that a single term can occur 
with more than one degree for its coefficient; in particular, the term (2 — g) occurs with both 
coefficient c^^^ and coefficient c^^^. Given that this phenomenon occurs, the sum (4.23) is 
not, in general, identical to the sum of distinct A's but instead is an upper bound on this 
sum. This is a type of behavior that does not occur for the Tutte or chromatic polynomials 
of strips of regular planar lattices that we have considered in previous work. For these strips, 
the sums of the corresponding numbers riT^Ly, d) and np{Ly, d) over d for a given value of Ly 
are equal to the respective numbers of distinct A's in the Tutte and chromatic polynomials. 
From our general calculation of the Tutte polynomial in [64] , we find that the nonzero terms 
are 



A 



Asg,3Ai = 9 - 5g + 7 
Asg,3,i,i — 2 — q 

10 - 6g + g^ ± (52 - 56g + 28g^ - 8g^ + q^f^ 

Asg,3,2,i = 2 - g 



for i = 2,3 



Asg,3,2,2 - 4 



A 



sg,3,3 



1 . 



(5.3.1) 
(5.3.2) 

(5.3.3) 

(5.3.4) 
(5.3.5) 
(5.3.6) 
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We label the terms in the flow polynomial for the Mobius strip as Xsq,Ly,d,±,j where d, ± 
means that the term has the coefficient ±0*^*^^ For the terms Xsq,Ly,d,±,j occurring in the ffow 
polynomial for the Mobius strip with width = 3 we find 



A 



sq,3,0,+,l 



A 



sq,3,0,l 



q-5q + 7 
2-q 



Kq,3,0,+,2 — Xsq,3,2,l 



sq,3,l,+,j 



, = A 



sq,3,l,j+l 



10-6g + g^±(52 - 56g + 28g^-8g^ + g^)^/2 for j = 1,2 (5.3.10) 



Kq,3,l,-,l — '^sq.S,!,! — 2 — 
As5,3,2,+,1 = Asg,3,3 = 1 . 



(5.3.7) 

(5.3.8) 
(5.3.9) 



(5.3.11) 

(5.3.12) 



Then 



and 



F{sq[3 X m, eye], q) = (A,,,3,o,i)"' + c^'^ J](A.,,3,ij)' 



F{sq[3 X m, Mb.],q) = (A,,,3A+,i)'" + (A.,,3,0,+,2)"^ - (A.,,3,0,- 



+cW((A 



sq,3,l,+,l) 



sq,3,l,+,2) 



-(A. 



,,3,l,-,l)'" ) + C^'^ 



(5.3.13) 



(5.3.14) 



We observe that F{sq[3 x m, ci/c.],g) and F{sq[3 x m, Mb.],q) always have the factors 
{q — l){q — 2). For m > 3 odd, F{sq[3 x m,cyc.],q) also has the factor {q — 3). The flow 
numbers for the cyclic 3 x m strip are therefore the same as for the cyclic 2 x m strip of the 
square lattice, given above in eq. (5.2.8). Further, we have 



We also find 



(j){sq[3 X m,Mb.]) = 3 . 

6 if m > 2 is even 



F{sq[3 X m, eye], q = 3) 



F{sq[3 X m,Mb.],q^ 3) 



if m > 1 is odd 
2 if m > 2 is even 

4 if m > 1 is odd 



(5.3.15) 

(5.3.16) 
(5.3.17) 
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Re(q) 



Figure 1: Singular locus B in the q plane for fl{sq,3 x oo, cyc./Mb.,q) for the 3 x oo strip of the square 

lattice with cyclic or Mobius (Mb.) boundary conditions. For comparison, zeros of the flow polynomial 
F{sq, 3 X Lx, eye, q) for = 30 (so that this polynomial has degree equal to c(G) = 61) are also shown. 
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-0.2 0.2 0.4 0.6 0.8 

Re(u) 



Figure 2: Singular locus B in the u plane for fl{sq,3x oo, eye. /Mb., q) for the 3 x oo strip of the square lattice 
with cyclic or Mobius boundary conditions. For comparison, zeros of the flow polynomial F{sq,3x L^, eye, q) 
for Lx = 30 are also shown. 
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The locus B for L^; — > cxo is shown in Fig. 1. Again, this locus is noncompact in the q 
plane. The locus is shown in the u plane, where it is compact, in Fig. 2. The locus separates 
the q plane into several regions. Three of these regions, Rj, j — 1, 2, 3, contain intervals of 
the real axis, which are > 3 for 2 < g < 3 for i?2, and g < 2 for i?3. Hence, 

qcf = 3 for {G} = sq, 3 x oo, cyc./Mb. (5.3.18) 

A general feature that we find is that, aside from the trivial case Ly = 1, for which i5 = 0, i3 
crosses the real axis at g = 2 and g = 3 for all of the widths Ly = 2, 3, 4 for which we have 
obtained exact general formulas for the flow polynomials (see below for the Ly = 4 case). For 
Ly — 3, the dominant terms in the above-mentioned three regions are (with an appropriate 
choice of branch cut for the square root) Asq,3,i,2 in Ri, Xsq,3,2,2 in R2, and Xsq,3,i,2 in R3, so 
that 

//(sg[3 X 00, cyc./Mb.], q) = {\sq,3,i,2f'^ for g G i?i . (5.3.19) 
In region i?2, Xsq,3,2,2 is dominant, so that 

|//(sg[3 X oo,cyc./M6.],g)| = |g-4|i/2 for g g i?2 • (5.3.20) 

In region R^, Xsq,3,i,2 is dominant, so that 

|/(sg[3xoo,q/c./M6.],g)|-|A,,,3,i,2r/' for g e i?3 . (5.3.21) 

There is also a pair of small complex-conjugate enclosed regions located to the upper and 
lower right of the point g = 3; in these regions, A^^ 3 2,1 is dominant. The part of B forming 
the boundary between region R2 and regions i?4, R*^ is the vertical line segment given by 
Re{q) = 3, — 1 < Im{q) < 1. Four curves on B meet at the intersection points g = ±i, and 
as one moves further upward above q = i, B consists of a boundary separating regions Ri 
and i?3. This boundary bifurcates at g ~ 2.6 -|- 1.8i, and as one moves to larger values of 
Im{q), one enters the region R^. Correspondingly, there is the complex-conjugate region R^. 
Expanding the equation for the degeneracy of magnitudes of leading terms A around u — 0, 
we have 

\l-5u + 7u^\^ \1 - 5u + Su^ + 0{u^)\ (5.3.22) 

Expressing this in terms of polar coordinates, with u — re*^, and expanding around r — 0, 
we get the equation, to leading order in r, 

cos 2^ = as r^O (5.3.23) 

The solution to this equation is 

9=^^i±Dl^ J = 0,1, 2, 3 (5.3.24) 
i.e., 9 — ±7r/4, ±37r/4, so that the curves approach u — with these angles. 
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5.4 Ly = 4 Cyclic and Mobius Strips of the Square Lattice 

For the flow polynomial of the cyclic 4 x m strip of the square lattice our general results yield 
71^(4,0) = 3, nF(4, 1) = 6, 71^(4,2) = 6, nF(4,3) = 3, 71^(4,4) = 1 with the total number 
-^F,4,A = 19. We have calculated the flow polynomial for this case. We flnd that the three 
^^^,4,0 J, j — 1, 2, 3, with coefficients c^^^ are the roots of the equation 

e - {q' - Sq^ + 24q - 26)^' 

-(q^ - Vlq^ + 59g^ - 14^ + 193g - 101)C + - 3)(? - 2)^ = . (5.4.1) 
For the A's with coefficient c*^^\ we have 

Af,4,i,i = -{q - 2)' (5.4.2) 

AF,4,i,2 = -(g-3)^ (5.4.3) 
The Ai;'^4_ij for j = 3, 4, 5, 6 are the roots of the equation 

- (g - 3)(g^ - 6g + 12)^^ - {q - 2){2q^ - 22q^ + 96?^ - 196? + 159)^^ 
-{q- 3)(g^ - 16g^ + 105g^ - 366g2 + 717g2 - 750g + 327)^ 

+ (g^ - 15g^ + 93g^ - 306?^ + 565g^ - 555g + 227) (g - 2)^ = . (5.4.4) 
For the A's with coefficient c^^\ 

Af,4,2J = \ {-{q - 3)' ± V(5-3)(g3-5g2 + llg-ll) j for j = 1, 2 . (5.4.5) 

The \F,A,2,j for j — 3, 4, 5, 6 are roots of the equation 

+ (2g2 - I3g + 22)i^ + (g - 2)(g3 - 13^^ + 51g - 66)^^ 

-{2q^ - 28?^ + 152?^ - 402g^ + 521? - 265)^ 

+ (g^ - llg^ + 43g2 - 70g + 41) (g - 2)^ = . (5.4.6) 

For the A's with coefficient c^^^ we flnd 

Af,4,3,i = g - 3 (5.4.7) 

Af,4,3,2 ^q-{?>-^/2) (5.4.8) 

Af,4,3,3 = g - (3 + \/2) (5.4.9) 
and, in accordance with eq. (4.11), 

Af,4,4 = -1 . (5.4.10) 
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Hence 



3 6 

F{sq[Axm,cyc.],q) = $^(A,,,4,o,,)™ + c^^^ $^(As,,4,i J'" 

3=1 j=l 

(5.4.11) 

For the Mobius strip of the square lattice with width = 4 we find 77,^(4,0,+) = 5, 
nj.(4, 0, -) = 4, 71^(4, 1, +) = 4, nj.(4, 1, -) = 3, 71^(4, 2, +) = 2, nj.(4, 2, -) = 1, with 



so that 



Asq,4,0,+j — Asq,4,0j , J — 1, 2, 3 

Asq,4,0,+,4 = Asq^4,2,l , Asq,4,0,+,5 = '^sq,4,2,2 
Asg,4,0,-,j = Kq,4,2,j+2 , J = 1, 2, 3, 4 
Asg,4,l,+j = ^sq,4,l,j+2 , i = 1, 2, 3, 4 

As<},4,l,-j = Asg^4^ij , j = 1,2 
Asq,4,l,-,3 = Asg^4^4 = —1 
Asq,4,2,+,1 = Asg,4,3,2 , Kq,4,2,+,2 = Kq,4,3,3 
Asq,4,2 ,-,1 = Kq,4,3,l 



F{sq[4 X m, Mb.], q) = J](A,,,4,o,+j)"' - Y.^^^^'^'^,- 
+c^'^[E(As,,4,i,+,)'"-E(A.,Ai,-.)™] 

2 



(5.4.12) 

(5.4.13) 
(5.4.14) 
(5.4.15) 

(5.4.16) 
(5.4.17) 
(5.4.18) 
(5.4.19) 



(5.4.20) 



The locus B for L^; — > cxd is shown in Fig. 3. This locus is again noncompact in the q 
plane, containing six curves that extend infinitely far away from q — 0. The locus separates 
the q plane into several regions. Three of these regions, Rj, j = 1,2,3, contain intervals of 
the real axis, which are g > 3 for 2 < g < 3 for i?2, and g < 2 for i?3. Hence, 

qcf = 3 for {G} ^sq, 4 x oo, cyc./Mb. (5.4.21) 

The locus B again crosses the real axis at g = 2 and q = 3. In region Ri the dominant term 
is the root of maximal magnitude of the cubic equation (5.4.1), which we denote Xsq,4,o,jmax, 
so that 

//(sg[4 X oo, cyc./Mb.], q) = {Xsq,4,o,jmaxy^^ for qeRi . (5.4.22) 
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Im(q)0 




Re(q) 



Figure 3: Singular locus B in the q plane for fl{sq, 4 x oo, eye. /Mb., q) for the 4 x oo strip of the square lattice 
with cyclic or Mobius boundary conditions. For comparison, zeros of the flow polynomial F{sq,4:X L^, eye, q) 
for Lx = 30 are shown. 
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In region R2, the dominant term is the root of maximal magnitude of the quartic equation 
(5.4.6), denoted \sq,4,2,jmax, so that 

\fl{sq[4: X 00, cyc./Mb.], q)\ = \Xsg,4,2,jmax\^^^ for q e R2 ■ (5.4.23) 

In region R^, the dominant term is the root of maximal magnitude of the other quartic 
equation (5.4.4), so that 

|//(sg[4 X 00, eye. 

/Mb.],q)\ = \Xs,,4,i,jmax\'^'' for g G i?3 • (5.4.24) 

In addition to the regions Rj, j = 1, 2, 3 that contain intervals of the real axis, there are also 
four complex-conjugate pairs of regions away from the real axis, Rj, R*, j = 4, 5, 6, 7. These 
can be identified in Fig. 3 as follows: i?4 is a small "bubble" region containing the point 
q = 3.1 + 0.6i; Rq contains the point g = 2 + 3i and extends to complex infinity; Ry contains 
the point q — Si and extends to complex infinity; and R5 is a very small region that contains 
the point g = 2.4 + 1.8i 

5.5 On Qcf for Cyclic/Mobius Square-Lattice Strips 

Aside from the trivial case Ly — 1, all of the cyclic/Mobius strips of the square lattice that 
we have studied have yielded qcf = 3. This is the same value as for the infinite square 
lattice (see eq. (17.1). We are led to conjecture that qcf = 3 for cyclic/Mobius strips of 
the square lattice for all Ly > 2. Our present finding and conjecture are related to our 
previous result that qc — 3 for Bw for self-dual cyclic strips of the square lattice [57, 68] for 
all widths considered and our consequent conjecture that this property qc — 3 holds for all 
cyclic self-dual strips of the square lattice. To see this connection, we recall that the dual 
graph of the cyclic strip of the square lattice with Lj^ > 2 is a cyclic strip of this lattice 
with width Ly — 1 augmented by edges joining all of the vertices on the upper and lower 
sides of the strip to two respective external vertices. The latter strip is rather similar to the 
cychc self-dual strips studied in [57, 68] , differing only by the edges joining the lower vertices 
to the second external point. To the extent that these lower edges do not modify qc, one 
would then expect that qc = 3 for the family of dual graphs, which is the equivalent to the 
conjecture that we make here for qcf. 

Our finding that, aside from the trivial case Ly = 1, B for the Ly x 00 cyclic/Mobius 
strips of the square lattice crosses the real axis at q — 2 for Ly — 2,3,4 suggests another 
conjecture, namely that this property holds for all cyclic or Mobius strips with width Ly > 2. 
It is interesting to relate this to our earlier finding that for this class of cyclic/Mobius strip 
graphs of the square lattice, Bw crosses the real axis at g = and g = 2 for all widths 
1 < Ly < 5 considered [32, 43, 50, 61]. One sees that for the widths Ly > 2 considered, 
Bfi and Bw both cross the real axis at g = 2, while Bw, but not Bfi, crosses at g = and 
Bfi, but not Bw, crosses at g = 3. For the self-dual strips considered in [57, 68], for which 
Bfi = Bw, there are crossings at g = 2, 3 but not g = 0. The actual behavior of Bfi = Bw in 
the Ly —>■ 00 limit, i.e. for the full square lattice, would presumably combine these features. 
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Tabic 2: Comparison of vahics of fl(sq{Ly),FBCy,q) with fl{sq,q) for 3 < (j < 10. For each vahic of q. 
the quantities in the upper hne are identified at the top and the quantities in the lower hne are the values 
of Rfi{sq{Ly), FBCy, q). The FBCj^ is symbolized as F in the table. 



Q 


fl{sq{l),F,q) 


fl{sq{2),F,q) 


fl{sq{3),F,q) 


fl[sq[A),F,q) 


!l\sq,q) 
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1.4142 


1 


1 


1.1740 


1.53960... 






0.6495 


0.6495 


0.7625 
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1.73205 


2 


1.7989 


1.9520 


2.3370(7) 
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3 


2.7248 


2.8827 


3.2510(10) 






0.9228 


rv 0001 

0.8381 


0.8867 




6 


2.236 


4 


3.6802 


3.8418 


4.2003(12) 






0.9523 


0.8762 


0.9146 




7 


2.449 


5 


4.6502 


4.8138 


5.1669(15) 






0.9677 


0.9000 


0.9317 




8 


2.646 


6 


5.6286 


5.7934 


6.1431(20) 






0.9767 


0.9163 


0.9431 




9 


2.828 


7 


6.6124 


6.7779 


7.1254(22) 






0.9824 


0.9280 


0.9512 




10 


3 


8 


7.5998 


7.7657 


8.1122(25) 






0.9862 


0.9368 


0.9573 





6 Approach of fl[sq^ Ly x oo, q) to fl{sq, q) 

It is of interest to use our exact calculations of the flows per face, fl{sq, Ly x 00, q) for the 
cyclic/Mobius strips of the square lattice to study how the values of this function approach 
the values for the infinite square lattice as the strip width Ly increases. This is done in Table 
2 for a range of q values. For this table, we define the ratio 

R„(s,(L,),FBC„,) = . (6.1) 

We use calculations of fl{sq, q) = W{sq, q) for the infinite square lattice obtained using 
Monte Carlo methods from [78, 79]. Again, we note that the case Lj^ = 1 is atypical since 
the fiow polynomial does not exhibit exponential growth with L^ and the number of faces is 
fixed at two, independent of L^ — m, in contrast to all of the strips with Ly > 2. (Hence, 
for Ly — 1 we do not list the Rfi values.) For > 2 we find that fl approaches the infinite- 
width value from below and this approach is not, in general, monotonic. For a given value 
of Ly > 2, the values of // are closer to those for the infinite lattice for larger q. In contrast, 
for cychc/Mobius strips and for a given value oi Ly, W calculated for the infinite strip of 
width Ly approaches the infinite- width value from above rather than from below [79]. 
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7 Strips of the Square Lattice with Torus and Klein 
Bottle Boundary Conditions 



It is of interest to study flows on lattice strip graphs that have doubly periodic boundary 
conditions. We carry out this study in the present section for the square lattice strips, 
including both the case of torus and Klein bottle (Kb.) topologies. 

7.1 Ly = 2 Strips of the Square Lattice with Torus and Klein Bottle 
B.C. 

This family involves double edges on each transverse slice of the strip. Although the chro- 
matic polynomial is insensitive to the presence of multiple edges, the flow polynomial is 
sensitive to this property, as is obvious since in general multiple edges allow more flows 
satisfying the conservation conditions at each vertex. Using our results in [64] we have 

A.,t,2,o,i = ^4 = g'-3g + 3 (7.1.1) 

where 

(where the second equahty holds for m > 2) , 

A.,t,2,i,i = g'-4g + 5 (7.1.3) 
A.,t,2,2 = 1 (7.1.4) 

and 

= (g2 - 3g + 3)™ + (g - l){q^ -Aq + 5)"^ + (g^ - 3g + 1) 

(7.1.5) 

F{sq[2xL,^m,Kh],q) = {Kt,2fl,iT + c^''\Kt,2,i,iT - I 

= (g2-3g + 3)"^ + (g-l)(g2-4g + 5)'"-l . 

(7.1.6) 

We have C F,sq,Ly=2,tor = (q* — 1)^ for the 2 x m torus strip and Cp^sq,Ly=2,Kh. = — 1 for the 
2x171 Klein bottle strip of the square lattice. Both flow polynomials have {q — 1) as a factor, 
so that the flow number = 2 for these families of graphs. 
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The locus B is noncompact in the q plane and separates this plane into four regions. This 
locus is identical to the locus for Bw in Fig. 5(a) of [38]. The region Ri includes the real 
segment q > 2, while R2 contain the segment q < 2. For this family, 

qcf = 2 for sq, 2 x 00, torus or Kb. (7-1-7) 

This demonstrates that for a given type of lattice strip and for a given strip width Ly, the 
imposition of different transverse boundary conditions leads, in general, to a different locus 
B and, for the case of periodic longitudinal boundary conditions where B crosses the real 
axis, can lead to different values of qcf. 

Above and below qcf there are two complex-conjugate phases that extend up and down 
to triple points, from which c.c. curves extend outward to complex infinity. In Ri, Xsqt,2,o,i 
is dominant, so 

fl{sq[2x 00, torus/ Kb.], q) = {q'^-3q + 3Y/^ for q e Ri . (7.1.8) 

where the notation "torus/Kb." means that the result holds for the strip with either torus 
or Klein bottle boundary conditions. In region R2, Xsqt,2,i,i is dominant, so 

\fl(sq[2x 00, torus /Kb.], q)\^\q^ -Aq + 5\^/^ for q E R2 . (7.1.9) 

In regions R^ and R^, Xsqt, 2,2 is dominant, so 

\fl{sq[2xoo,t(yrus/Kb],q)\^l for q e R3, R^ . (7.1.10) 



7.2 Ly = 3 Strips of the Square Lattice with Torus and Klein Bottle 
B.C. 

Here, for the Ly — 3 strips of the square lattice with torus and Klein bottle (Kb.) boundary 
conditions, which we denote generically as sqSt and sq3kb, we find 

NF,sq3t,X = 8 (7.2.1) 

NF,sq3kb,X = 6 . (7.2.2) 

These are reductions from the numbers of different terms in the Tutte polynomials for these 
strips, which are [64] Nx^sqst.x = 20 and N'r^sq3kb,x = 12. These may be compared with the 
numbers for the chromatic polynomials for these strips [46] Np^sq3t,x = 8 and Np^sq3kb,x — 5. 
We have 

8 

F{sq[3 X m, torus], q) = Csq3t,j{Kq3t,j)"' (7.2.3) 

6 

F{sq[3 X m, Kb.],q) = J] Csq3kb,j{Xsq3kb,jr (7-2.4) 
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where, in order of increasing degrees of the coefficients (see below) 



A 



sqStJ 





= Q 


-2 


Kq3t,5 


= Q 


- 1 


^sq3t,6 


= Q 


-4 


^sq3t,7 


= Q 


-5 


^sq3t,8 




-1 . 



A.,3t,i = - 6?" + Uq - 13 (7.2.5) 

^[-18+16g-6g2 + g3±y^j for J = 2, 3 (7.2.6) 

Rsq3t = 256 - 440? + 376g^ - 196g=^ + 64g^ - 12?^ + (7.2.7) 

(7.2.8) 
(7.2.9) 
(7.2.10) 
(7.2.11) 
(7.2.12) 

In contrast to the situation with cychc and Mobius strips, there is not a 1-1 correspondence 
between terms A and the coefficients that the flow polynomial inherits as a special case of 
the Tutte polynomial; here, one of the terms appears with two different Tutte coefficients. 
Specifically, the term Xsq3t,4: appears with both the coefficient 2{q — 1) and q{q — 3), so that 
its net coefficient is the sum, {q + l){q — 2). For the coefficients we thus have 

Csq3t,l = 1 (7.2.13) 
Csq3t,2 = Csq3t,3 = ? " 1 (7.2.14) 

Csq3tA = {.q + l){.q-2) (7.2.15) 
Csq3t,b = ^C,g3t,6 = ^(? - !)(? - 2) (7.2.16) 

c.,3t,7 = ^g(g - 3) (7.2.17) 

Csq3t,s = - 6g2 + 8? - 1 . (7.2.18) 

The relation of these with the c'^'^^'s was discussed earlier [52, 61]. The sum of coefficients is 
CF,sq,Ly=3,tor = — 1)^ for the 3 X m torus strip. 
For the Klein bottle strip we have 

Kq3kb,j = Xsq3t,j , J = 1; 2, 3 (7.2.19) 

Kq3kb,4: — Kq3t,5 — Q — ^ (7.2.20) 
Kq3kb,5 — Kq3t,7 — q — 5 (7.2.21) 
Kq3kb,6 — Kq3t,8 — "1 ■ (7.2.22) 

Thus the flow polynomial for the Ly — 3 strip of the square lattice with torus boundary 
conditions has two terms, q — 2 and q — 4:, that are absent from the flow polynomial of this 
lattice strip with Klein bottle boundary conditions. Evidently, the set of terms for the strip 

30 



with Klein bottle boundary conditions is a subset of the set of terms for the strip with torus 
boundary conditions. 

The corresponding coefficients are 

Csq3kb,l = Csq3t,l = 1 (7.2.23) 



Csq3kb,j = Csqst,j = ? " 1 , J = 2, 3 (7.2.24) 
1 

2I 



1 

Csq3kbA = Csg3t,5 = 7^(? " 1)(? " 2) (7.2.25) 



Csq3kb,5 = Csq3t,7 = ^^(9 " 3) (7.2.26) 

Csq3kb,e = ~iq-l) ■ (7.2.27) 

The sum of these coefficients is CF,sq,Ly=3,Kb. = (q — 1)^- The flow numbers of the sq, 3 x 
strips with torus and Klein bottle (kb) boundary conditions are 

(/)(sg[3 X L^, torus/ Kb.]) = 2 . (7.2.28) 

The locus B is the same for the torus and Klein bottle strips of a given lattice. We have 
proved this type of equality for Bw in [52] for chromatic polynomials of torus and Klein 
bottle strips. The proof for the chromatic polynomials uses the construction of a family of 
lattice strip graphs which are identical to the torus and Klein bottle strips for even and odd 
length L^. The locus B is the accumulation set of the zeros in the limit L^, oo, and this 
limit can be taken over either even or odd values of L^, which proves the equality of the loci 
B for the respective Lj; — > oo hmits of these two types of strips. The same method of proof 
applies here for the flow polynomials. Note that this result requires that all of the dominant 
terms arc the same for the strips with torus and Klein bottle boundary conditions. This 
condition is satisfied. 

The locus B for the ^ oo limit of the Ly = 3 square strip with torus or Klein bottle 
boundary conditions is shown in the q and u planes in Figs. 4 and 5. As is evident, B 
is noncompact in the q plane and separates this plane into various regions. Three of these 
regions, Rj, j = 1, 2, 3, contain intervals of the real axis, which are q > qcf for Ri, 2 < q < qcf 
for i?2, and g < 2 for i?3, where 

gc/ = 2.6120932... for sg, 3 x oo, torus/Kb. (7.2.29) 

This is the larger of the two real roots of the equation 

2g^ - 19g^ + 71q^ - U2q + 132 = (7.2.30) 

resulting from the condition of degeneracy in magnitude of the leading terms lA^^stjl = 
I'^sg3t,2|- The fact that qcf increases from the value 2 for Ly = 2 to approximately 2.61 
for Lj^ = 3 is in accord with the property that as Ly — > oo, qcf — > 3 (see eq. (17.1) 
below). The behavior exhibited by the strips of the square lattice with torus/Klein bottle 
boundary conditions is thus that qcf progressively approaches the asymptotic value qcf — 3 
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Figure 4: Singular locus B in the q plane for fl{sq,3 x oo, torus/ Kb., q) for the 3 x oo strip of the 
square lattice with torus or Klein bottle boundary conditions. For comparison, zeros of the flow polynomial 
F{sq, 3 X Lx, torus, q) for = 30 are also shown. 
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Figure 5: Singular locus B in the u plane for fl(sq,3 x oo, torus/Kb. ,q) for the 3 x oo strip of the 
square lattice with torus or Klein bottle boundary conditions. For comparison, zeros of the flow polynomial 
F{sq, 3 X Lx, torus, q) for = 30 are also shown. 
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as Ly increases rather than being equal to this value for each > 2 as we find for the 
cychc/Mobius strips of the square lattice. 

The dominant term in region Ri is Xsq3t,2 so that 

fl{sq[3 X oo,torus/Kb.],q) = (A,,3t,2)'/' for q e Ri . (7.2.31) 

This is interesting since, unlike the flow polynomials for lattice strips that we have studied 
before, here the dominant term in region Ri does not have coefficient 1. We recall that 
the chromatic polynomials for the recursive strips that we have considered always have 
the property that the dominant term in region Ri (including the positive real q axis) has 
coefficient 1 [47]. In region R2, Xsqstj = q — 5 is dominant, so that 

\fl{sq[3 X 00, torus/ Kb], q)\ = |g - 5p/^ for q e R2 ■ (7.2.32) 

In region R^ Kq3t,3 is dominant, so 

\fl{sq[3x 00, torus/ Kb.], q)\^\X,,st,3\^/^ for q e R^ . (7.2.33) 

The boundary on B separating regions Ri and R2 is the vertical line segment in the q plane 
given by Re{q) = qcf as given above in eq. (7.2.29) and —1.04 < Im{q) ^ 1.04. This 
line segment ends at the complex- conjugate set of triple points at g = qcf ± 1.04i. Going 
vertically upward and downward away from the real axis along this direction, one passes into 
two complex-conjugate regions, i?4 and RI in which \sq3t,i is dominant, so that 

\fl{sq[3x cx), torus /Kb.], q)\ = \q^ -Gq"^ + 14:q-13\^/^ for e R4, Rl . (7.2.34) 

At the point u = 0, four curves on B intersect. Let us express the terms Xsq3t,j in terms of the 
variable u = 1/q and defining the scaled terms Xsq3t,j,u = u^Xsq3t,j. With appropriate choices 
of branch cuts for the square root, the equation of degeneracy of magnitude of leading terms 
in the vicinity oi u — 0, when expanded in a small u series, is 

]l-6u + 15w^ + 0{u^) \^\l-6u + Uu'^ + 0{u^) ] . (7.2.35) 

Using polar coordinates u — re*^, this reduces to eq. (5.3.23) as r — > 0, which shows that 
the curves approach u — with the angles 9 — ±7r/4, ±37r/4. 

8 Self-Dual Strips of the Square Lattice 

In [57, 68] we calculated chromatic polynomials and Tutte polynomials for families of planar 
self-dual strips of the square lattice with free and cyclic longitudinal boundary conditions, 
denoted, respectively, Go^LyXL^, eye.) and Go^LyXL^, free). Applying the relation (3.1.2), 
we infer, in particular, that 

F{Gn,q)^q-'P{Gu,q) (8.1) 
where G d refers to any of these graphs. Hence, 

Qcf — Qc — ^ for Gi) famihes . (8.2) 

Similarly, the loci B are the same for the L^, — > 00 limits of the flow and chromatic polyno- 
mials for these families of graphs. A particular property is that, in contrast to the loci B for 
the other families studied here, these are compact in the q plane. 
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9 Strips of the Square Lattice with Free Boundary 
Conditions 



It is of interest to compare the flow polynomials for strips with periodic or twisted peri- 
odic longitudinal boundary conditions to those with free boundary conditions. These flow 
polynomials may be calculated directly by iterative application of the contraction-deletion 
property or via the chromatic polynomials of the dual graphs, or as special cases of Tutte 
polynomials. We list some relevant results. Clearly for a tree graph, the flow polynomial 
vanishes. For strips with free longitudinal boundary conditions, we shall use a different 
labelling convention, viz., m = — I, than the convention m — used for strips with 
periodic longitudinal boundary conditions. For m > 1 we have 

F{sq[Ly = 2, m + 1, free], q) = {q - l){q - 2)— ^ (9.1) 

F{sq[Ly = 3, = m + 1, free], q) = {q - l)(g - 2){q^ - 5q + 7)"^-^ . (9.2) 

Thus, NF^sq,Ly,\ = 1 for Ly = 2, 3. For both of these families of strip graphs, the continuous 
accumulation set of zeros is the empty set since the zeros are discrete. 

It is convenient to use a generating function to give the results for the cases Ly > 4. This 
is similar to the generating functions that were utilized for chromatic polynomials in [34, 36] 
and our subsequent works. For the strip of type G, length = m + 1, we write 

oo 

r(G,?,^) = ^F(G^,g)z- (9.3) 

where 



m=0 



with 



H{G,q,z) = Y,MA<D^' (9-5) 

and 

V{G,q,z) = l + J2bG,Mz' (9.6) 

Af{sq[Ly = 4, free], q, z) = {q- l){q - 2)z[{q -2)-{q-l){q- ?,)z] (9.7) 



V{sq[Ly = 4, free],q, z)^l-{q^ - %q^ + 24g - 26)z 

-(g^ - 12g^ + 59g^ - 149?^ + 193g - 101)^^ + (g - 2f{q - 2,)z^ . (9.8) 

Thus, NF^sq,Ly=ijree,\ = 3. The locus B for the — > oo limit is compact and consists of 
complex-conjugate pairs of arcs together with a self-conjugate arc that crosses the real axis 
at g ~ 2.6252. This is similar to what we found for the loci Bw for strips with free boundary 
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conditions in earlier work [34, 36]. Although B crosses the real axis for this strip, it is not 
guaranteed to cross the real axis when one uses free longitudinal boundary conditions, as 
will be illustrated by an explicit example below. 

We have also calculated the generating function for F{sq[Ly — 5, L^, free\,q) and find 

that NF,sq,Ly=5,free,X = 4, with 

^sq,5Jree,0 = (9.9) 
^,,5,/.ee,l=(g-l)(g-2)3 (g^Q) 

Asq,5,free,2 = - - 2)(g^ - 8?^ + 22q^ - 20q^ - 5g + 9) (9.11) 

Asq,5,free,3 = (?-!)(?- 2) {2q' - 29g^ + 178q' - 598g^ + 1186g3 - 1387?^ 

+884g - 237) (9.12) 

Asq,5,free,4 = -(? -!)(?- 2)^(g^ - Uq' + 79g^ - 229g3 + 359q^ - 288g + 91) (9.13) 
bsq,5,free,l = + 10?^ - Uq^ + 97q - 88 (9.14) 
bsq,5,free,2 = + 19?^ - 156?^ + 726?^ - 2085g^ + 3711?^ - 3790g + 1708 (9.15) 

bsq,5,free,3 = (? " 3) {2q^ - 39g^ + 337g^ - 1685?^ + 5335g^ - 10959g^ + 14264?^ 

-10753g + 3595) (9.16) 

bsq,5,free,4 = -{q " ^fiq^ - 23q^ + 235^^ - 1401g^ + 5376g^ - 13785g^ + 23647^^ 

-26198?^ + 17033g - 4964) . (9.17) 

The locus B is comprised of arcs and is similar to, although more complicated than, that for 
the Ly = 4 strip. 

10 Strips of the Square Lattice with Cyhndrical Bound- 
ary Conditions 

For strips of the square lattice with cylindrical boundary conditions we calculate, for m > 1, 
F{sq[Ly^2,L,^m+l,cyl.],q) = (g - l)(g - 2)2(L>4)— ^ 

= (q-l)(q-2nq'-3q + 3r-' (10.1) 

F{sq[Ly = 3, = m + 1, cyl.], q)^{q- l){q -2){q- 3f{q' - 6q^ + Uq - 13)'""^ . (10.2) 

For both of these families of strip graphs, the continuous accumulation set of zeros is empty 
since the zeros are discrete. 
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The coefficients of the generating function for F{sq[Ly = 4, L^, cyl.], q) are 

Asq,4,cyl,0 = ? " 1 (10-3) 

Asq,4,cyi,i = -{q - l)(3g' - 18?' + 35g - 18) (10.4) 

Asg,4,cyi,2 = {q- lf{q - 3)(g=^ - + 13q - 9) (10.5) 

bsgA,cyi,i = -Q^ + - 29g2 + 55g - 46 (10.6) 

bsq,4,cyi,2 = q^ - I2q^ + 61g^ - 169g^ + 269g' - 231g + 85 . (10.7) 
Thus, NF,sq,Ly=4,cyi,x = 2, and the terms are 

1 



^F,sq,Ly=4,cyl,j 



2 



- 8q^ + 29?' - 55g + 46 ± ^/R^c (10.8) 
where 

Rsq4c = - 16?^ + 118?^ - 526g^ + 1569g^ - 3250?^ + 4617g' - 4136g + 1776 . (10.9) 

The denominator I) of the generating function is the same as that for the generating function 
for chromatic polynomials of the Ax L^^ strip of the square lattice with cylindrical boundary 
conditions, and consequently, in the L,,. oo limit, the locus B is also the same as the 
locus Bw for the 4 x oo square-lattice strip with cylindrical boundary conditions, shown as 
Fig. 3(a) in [36]. This locus is compact and consists of a complex-conjugate pair of arcs 
together with a self-conjugate arc that crosses the real axis (at q ~ 2.30) and a very short 
line segment emanating outward from this crossing point on the real axis. The value of qcf, 
which is given by the right-hand end of the very short line segment, occurs at g ~ 2.35. The 
locus does not separate the q plane into different regions. 

The coefficients of the generating function for F{sq[Ly = 5, L^, cyl.], q) are 

Asq,5,ayl,0 = ? " 1 (10.10) 

Asq,5,cyi,i = -(g - l)(4g^ - 35g3 + 115g' - 157g + 62) (10.11) 

Asq,5,cyi,2 = "(? " !)'(?' " 7q' - + " 458?' + 645g - 348) (10.12) 

bsqAcyLi = -q^ + lOg' - 46g^ + 124g2 - 198g + 148 (10.13) 

bsq,^,cyi,2 = - 19?^ + 159g^ - 767g^ + 2339g^ - 4627g^ + 5800g' - 4212g + 1362 . (10.14) 

Thus, Nsg,Ly=5., ,cyi.,\ — 2. As was the case with Ly — 4, the denominator I? of the generating 
function is the same as that for the generating function for chromatic polynomials of the 
5 X strip of the square lattice with cylindrical boundary conditions, and consequently, in 
the Lx ^ oo limit, the locus B is also the same as the locus Bw for the 5 x oo square- lattice 
strip with cylindrical boundary conditions, shown as Fig. 2 of [50]. As before, B is compact 
and consists of five arcs, four of which form two complex-conjugate pairs and one of which is 
self-conjugate. The self-conjugate arc crosses the real axis at qcf — 2.69168. The endpoints 
of the arcs occur at the ten zeros of the polynomial i)lq,b,cyi,i " ^bsg,5,cyi,2, which square roots 
that occur in Xsq,Ly=5,cyi,j have branch point singularities. The arcs comprising B do not 
separate the q plane into different regions. 
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11 Ly = 2,3 Cyclic Strip of the Honeycomb Lattice 

The Ly = 2 cyclic or Mobius strip of the honeycomb lattice can be constructed by starting 
with the Ly = 2 cyclic or Mobius strip of the square lattice and inserting degree-2 vertices 
on each horizontal (longitudinal) edge. Because the insertion of degree-2 vertices does not 
affect the flow polynomial, it follows that, for BC^ — FBC^, PBC^, TPBC^ and Ly — 2, we 
have 

F{hc[2 X L^, FBCy, BC^i q) = F{sq[2 x L^, FBCy, BC^^q) (11.1) 

and 

fl{hc[2 X oo, cyc./Mb], q) = fl{sq[2 x oo, cyc./Mb],q) (11.2) 

where the expressions for F{sq, 2xLx, FBCy, BC^, q) with the various longitudinal boundary 
conditions were given above. 
For = 3 we calculate 

3 

F{hc[3xm,cyc.],q) = (Afe,,3,o,i)'" + c^^' 5^(Afee,3,ij)" 

+ C^'^J2(>^Hc,S,2,r + C^'^ (11.3) 



where 



A^c,3Ai = (?-3)' (11.4) 
Xhc,3,2,j = l[7 - 2q± ^13 - Aq]^ for j = l,2 (11.5) 

A;,c,3,3 = 1 (11.6) 
and the Xhc,3,i,j, i = 1, 2, 3 are the roots of the equation 

e-iq- 3)(g - -{q- Sf{2q - 5)^ - {q - 2)\q - 3)^ = . (11.7) 

Hence, 

0(/ic[3 X m,c?/c.]) = 4 . (11.8) 

The locus B for the — oo limit of the cyclic or Mobius Ly = 3 strips of the honeycomb 
lattice is shown in the q plane in Fig. 6 and in the u plane in Fig. 7. We have 

5 + \/5 

qcf = — ^ = 3.61803... for he, 3 x oo, cyc./Mb. (11.9) 

This is already within 10 % of the asymptotic value q^f — 4 for the 2D honeycomb lattice 
(see eq. (17.2) below). The locus B is again noncompact in the q plane, passing through the 
origin of the u plane. The locus separates the q plane into several regions, which contain 
intervals of the real axis: j = 1, 2, 3, 4, which include the respective intervals (1) q > qcf, 



38 



Re(q) 



Figure 6: Singular locus B in the q plane for fl{hc, 3 x oo, eye./ Mb, q) for the 3 x oo strip of the honeycomb 
lattice with cyclic or Mobius boundary conditions. For comparison, zeros of the flow polynomial F{hc, 3 x 
Lx, eye, q) for = 30 are also shown. 
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Figure 7: Singular locus B in the u plane for fl{hc, 3 x oo, eye. /Mb, q) for the 3 x oo strip of the honeycomb 
lattice with cyclic or Mobius boundary conditions. For comparison, zeros of the flow polynomial F{he,3 x 
Lx, eye, q) for = 30, expressed in terms of u, are also shown. 
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(2) 3 < g < Qcf, (3) 2 < g < 3, and (4) q < 2. In region Ri, the dominant term is the root 
of the cubic equation (11.7) with maximal magnitude, which we denote as Xhc,3,i,i, so that 



fl{hc[3 X oo, cyc./Mb],q) = (Afee,3,i,i)'^' for q e Ri . 
In region R2, Xhc,3,3 = 1 is dominant, so 

\fl{hc[3 X 00, cyc./Mb], q) \ = I for q e R2 . 
In region R^, Xhc,3,2,i is dominant, so 



\fl{hc[3 X 00, cyc./Mb], q)\ 



^(7-2g+ Vl3-4g 



1/2 



for q e R3 . 



;ii.io) 



(11.11) 



:ii.i2) 



In region R4 the maximal root of the cubic equation (11.7) is dominant, with a result 
analogous to (11.10) for //. In regions R5 and R^ the dominant A is \hc,3,o,i, so that 

\fl{hc[3x 00, cyc./Mb], q)\^\q-3\ for q e R5, R^ . (11.13) 

The curves cross the origin of the u plane with the angles ±7r/4 and 37r/4. 

12 Ly = 2 Cyclic Strip of the Triangular Lattice 

Next, we consider strips of the triangular lattice. For the Ly — 2 cyclic strip we calculate 

AtH,2,o,i = - 2)' (12.1) 

1 



Xtri,2,ij = - 6 - 4g + g2 ± (g _ 2)^8 -4q + q^ 



,3 2 

where j = 1, 2 for the ± sign before the square root, and 



Xtri 



tri,2,2 



1 . 



Then 



(2) 



(12.2) 



;i2.3) 



(12.4) 



F{tri[2 X m, eye], q) ^ {q - 2f^ + c^^) [(AtH,2,i,i)" + (AtH,2,i,i)"J + c 

Thus, NF^tri,Ly=2,cyc,\ — 4, and for the cyclic strip of the triangular lattice of width Ly — 2, 

nF{tri,2,0) ^ 1, ^^(iri, 2, 1) = 2, npitri, 2,2) ^ 1 (12.5) 

with npitri, 2,d) = Q for d > 3. For this strip the sum of coefficients is 

CF,tri,Ly=2 = q{q - 1) for tri, 2 x L^, eye. (12.6) 

Using our exact solution, we observe that F{tri[2 x ■m,cyc.],q) has only the common 
factor (g — 1) and hence 

(f){tri[2 X m, eye.]) = 2 . (12.7) 



41 




Figure 8: Singular locus B in the q plane for fl{tri, 2 x oo, eye. /Mb, q) for the 2 x oo strip of the triangular 

lattice with cyclic or Mobius boundary conditions. For comparison, zeros of the flow polynomial F{tri, 2 x 
Lx,cyc.,q) for = 30 (so that this polynomial has degree 61) are also shown. 
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Figure 9: Singular locus B in the u plane for fl{tri, 2 x oo, eye. /Mb, q) for the 2 x oo strip of the triangular 

lattice with cyclic or Mobius boundary conditions. For comparison, zeros of the flow polynomial F{tri,2 x 
Lx, eye, q) for = 30, expressed in terms of u, are also shown. 
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Further, we find that 

(l){tri[2 X m, Mb.]) = 3 . (12.8) 

The locus B for the cyclic and Mobius Ly = 2 strips of the triangular lattice is shown 
in the q plane in Fig. 8 and in the u plane in Fig. 9. This locus is noncompact in the q 
plane and divides this plane into four regions, Ri, R2, and the complex-conjugate regions 
i?3 and i?3. Regions Ri and R2 contain the respective real intervals q > 2 and q < 2 and are 
contiguous along a vertical line segment extending from q — 2 + ito2 — i. Triple points on 
B occur at these endpoints q = 2±i. Evidently, 

qcf = 2 for tri, 2 x 00, eye./ Mb. . (12.9) 

Along the vertical line segment from 2 + i to 2 — i, |Atri,2,i,i| = |Atr-i,2,i,2| = ^^^,2,2 = 1- 
Regions Rz and i?3 extend upward and downward, respectively, from the complex-conjugate 
triple points aX, q — 2-\-i and q — 2 — i. In region \tri,2,\,i is dominant so that 

fl(tri[2 X 00, q/c], g) = f{tri[2 x 00, Mb.], q) = (AtH,2,i,i)'^' for q e Ri . (12.10) 

In region R2, Xtri,2,i,2 is dominant, so that 

\fl{tri[2x 00, eye./ Mb.], q)\ = \Xtri,2,i,2\'^^ for q e R2 . (12.11) 
In regions Rs and R^, Xtri,2,o,i is dominant, so that 

\fl{tri[2x 00, eye./Mb.],q)\ = \q-2\ for q E R3, RI ■ (12.12) 

Some special values include fl{2) = 1, |//(1)| = [(3 + V5)/2]^/\ and |//(0)| = [3 + 2V2]^/\ 
To show that B passes through -u = 0, we calculate the corresponding scaled terms 
Xtri,2,d,j,u = Xtri,2,d,j / (f with u = 1/q and see that the degeneracy condition of dominant 
Xtri,2,d,j,u^ has a solution at w = 0. The dominant Xtri,24j,u^ in the vicinity oiu — Q are 



A 



Atom = (1 - 2«)' (12.13) 

1 



tri,2,l,l,u 



1 - 4m + + (1 - 2u)\/l - 4m + 8^2 . (12.14) 
For ^ 0, the latter term has the Taylor series expansion 

At„,2,i,i,n = 1 - 4m + - M^ + O(w') . (12.15) 

Thus, in the same way as before, introducing polar coordinates and expanding the equation 
of the degeneracy of magnitudes of leading terms in the vicinity of m = yields the condition 
cos 29 — 0, thereby showing that four branches of B approach the origin of the u plane at 
the angles given by (5.3.24), i.e., 9 = ±7r/4 and 9 = ±37r/4. 
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13 Ly = 3 Cyclic Strip of the Triangular Lattice 

For the Ly — 3 cyclic strip of the triangular lattice we find the following results: 



A. 



tri,3,0,j — 



- 7q^ + 21q^ - 33q + 23± 



(g^ - 14g^ + 91q^ - 360?^ + 949?^ - 1708?=^ + 2047g^ - 1486? + 497) 



1/2 



J = 1,2 

(13.1) 
(13.2) 



= - 4g + 5 . 

The Xtri,3,i,j for j = 2, 3, 4 are the roots of the equation 

- (g^ - 7q^ + 22g^ - 37g + 30)^^ 
+{q^ - llq^ + h2q^ - l3Aq^ + 200?^ - 168g + 69)^ - {2q^ - Ihq^ + 42^^ - 51g + 24) = 

(13.3) 

(13.4) 



Atri,3,2,i = - 3g + 3 



A 



tri,3,2,j 



q'^ -5q + 9±{q^ - lOq^ + 43q'^ - 90q + 73) 



1/2 



, J = 2,3 



Atri,3,3 — 1 ■ 



(13.5) 
(13.6) 

Hence, NF,tri,Ly=z,cyc.,\ = 10, npitri, 3, 0) = 2, npitri, 3, 1) = 4, npitri, 3, 2) = 3, npitri, 3, 3) = 
1, and 

2 4 

F{tri[3 X m, eye], g) = ^(At,,,3,oj)"' + c^'^ ^(AtH,3,ij)"+ 

c(^)E(W2.i)"^ + ^^'^- (13-7) 

Using our exact result (13.8) we observe that F(tri[3 x m, eye], q') has only the common 
factor (g — 1), so that 

(t){tri[3 X m, eye.]) = 2 . (13.8) 

This can also be derived as a corollary of the basic theorem discussed in section (3.2), that 
a bridgeless graph admits a 2-flow if and only if all of its vertex degrees are even; here all of 
the vertices of the tri[3 x m, eye] are even. 

For the cyclic strip of the triangular lattice of width Lj, = 3, the sum of coefficients is 



CF,tri,Ly=3 ^q{q-^f for tri, 3x L^, eye. 



;i3.9) 



The locus B for the cyclic and Mobius Ly — 3 strips of the triangular lattice is shown in 
the q plane in Fig. 10. This locus is noncompact in the q plane, containing eight curves that 
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Figure 10: Singular locus B in the q plane for fl{tri,3 x oo, eye. /Mb, q) for the 3 x oo strip of the 

triangular lattice with cyclic or Mobius boundary conditions. For comparison, zeros of the flow polynomial 
F{tri, 3 X Lx, eye, q) for = 20 (so that this polynomial has degree 81) are also shown. 
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extend infinitely far away from q = 0. The locus separates the q plane into several regions. 
Three of these regions, Rj, j = 1, 2, 3, contain interval of the real axis, which are q > qcf for 
Ri, and 2 < q < qcf for i?2, and ^ < 2 for R3, where 



qcf = 2.213548 for tri, 3 x 00, cyc./Mb. . 



;i3.10) 



This is within about 15 % of the asymptotic value g^/ — 2.618 for the 2D triangular lattice 
(see eq. (17.3) below). In regions Ri and R3, the dominant term is the root of maximal 
magnitude of the cubic equation (13.3), which we denote Xtri,3,i,jmax, so that 



fl{tri[3 X 00, cyc./Mb.], q) = (At^.s.ijmaa;)^^^ for q G Ri. 
In region R2, Kri,3,2,2 is dominant, so that 

\fl{tri[3 X 00, cyc./Mb.], q)\ = \Xtri,3,2,2\^^'^ for q & R2 ■ 

In region R^, 

\fl{tri[3 X 00, cyc./Mb.], q)\ = \Xtri,3,i,jmax\^^^ for q e -^3- 



(13.11) 



(13.12) 



(13.13) 



In addition to the regions Rj, j = 1,2,3 that contain intervals of the real axis, there are 
also three complex-conjugate pairs of regions away from the real axis, Rj, R*, j — 4,5,6. 
These can be identified in Fig. 10 as follows: R4 contains the point q = 4 + 3i and extends 
to complex infinity; R^ contains the point q = 1.5 + 3i and extends to complex infinity; and 
Rq contains the point q = — 1 + 3i and extends to complex infinity. 



14 Ly = 2 Strip of the Triangular Lattice with Toroidal 
and Klein Bottle Boundary Conditions 

For the Ly = 2 strips of the triangular lattice with torus and Klein bottle boundary condi- 
tions, we find that 



Xfri 



tri2t,j 



ll-19? + 15g^-6?3 + ?^ 



±(129 - 446g + 727g2 - 722^^ + 479g^ - 218g^ + 66g° - 12q' + g= 



J = l,2 
(14.1) 



A 



tri2t,j 



14 - 20q + 15g2 - 6q^ + 



±(212 - 600? + 852q' - 776g'^ + 493?* - 220q^ + 66q^ - 12q^ + q 



7 I „8\l/2 



i = 3,4 
(14.2) 
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and 

At„2t,5 = 2 . (14.3) 

The corresponding coefficients are 

CtH2tj = l i = l,2 (14.4) 

Ctri2t,j = q - i = 3,4 (14.5) 
and ^ 

ctri2t,5 = -^q{q - 3) . (14.6) 

It follows that 

F{tri[2 X m,torus],q) = {\tri2t,i)"' + {Kri2t,2)'^ 

+{Q - 1) [{Ki2t,3r + {Ki2tAr] + lq{q - 3)2"^ . (14.7) 

The locus B for the torus and Klein bottle Ly = 2 strips of the triangular lattice is shown 
in the q plane in Fig. 11 and in the u plane in Fig. 12. This locus is noncompact in the 
q plane, and divides this plane into various regions. Three of these regions, Rj, j — 1, 2, 3, 
contain intervals of the real axis, which are q > qcf ior Ri, 2 < q < qcf for R2, and < 2 for 
i?3, where 

gc/ = 2.307144568... for tri, 2 x 00, torus/ Kb. . (14.8) 
The dominant terms in regions Rj, j — 1,2,3 are Xtri2t,i, \ri2t,b and Xtri2t,i so that 

fl{tri[2xoo,torus/Kh],q) = {Xtri2t,if'^ for q e Ri (14.9) 
\fl{tri[2xoo,torus/Kh.],q)\^2^'^ for q e R2 (14.10) 

\fl{tri[2 X ^,torus/Kh],q)\ = \\tri2t,^\^'^ for g G i?3 • (14.11) 

There are also three complex-conjugate pairs of regions, Rj, R*, j = 4,5,6. These can be 
identified in Fig. 11 as follows: R4 is a "bubble" region centered at the point q — 1.8 + 1.5i; 
i?5 is a small "bubble" region centered at the point g = 1.2 + i; and Rq contains the point 
q — 2i and extends to complex infinity. The dominant terms in regions Rj, j — 4, 5, 6 are 

Xtri2t,l, Xtri2t,5 ^ud Xtri2t,l- 

To show that B passes through m = 0, we calculate the corresponding scaled terms 
Xtri2t,j,u = Xtri2t,j/q^ with u — 1/q and see that the degeneracy condition of dominant 
Xtri2t,j,u^ has a solution at w = 0. The dominant terms in the vicinity oi u — Q are Xtri2t,i,u 
and Xtri2t,z,u- For w — > 0, these terms have the Taylor series expansion 

Xtri2t^,u = 1 - 6ii + Ibu^ - 19w^ + 0{u^) (14.12) 

Xtri2t,z,u = 1 - 6ii + 15^/2 - 20^3 + Oiu") . (14.13) 
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Figure 11: Singular locus B in the q plane for fl{tri,2x oo, torus /Kb., q) for the 2xoo strip of the triangular 
lattice with torus or Klein bottle boundary conditions. For comparison, zeros of the flow polynomial F{tri,2x 
Lx, torus. ,q) for = 20 are also shown. 
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Figure 12: Singular locus B in the u plane for fl{tri,2x oo, torus /Kb., q) for the 2 x oo strip of the triangular 

lattice with torus or Klein bottle boundary conditions. For comparison, zeros of the flow polynomial F{tri,2x 
Lx, torus, q) for = 20, expressed in terms of u, are also shown. 
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Hence, introducing polar coordinates u = re*^, the condition of degeneracy of magnitudes in 
the vicinity oi u — yields the condition 

r^cos^(3-4cos2^) = as r^O. (14.14) 

This proves that four branches of B approach the origin of the u plane at the angles 

(2/c + l)7r ^ i<^<6 (14.15) 

6 

i.e., at ^ = ±7r/6, 9 = ±n/2 and 9 = ±5n/6. 

15 Strips of the Triangular Lattice with Free Boundary 
Conditions 

In this section, we give the flow polynomials of the triangular lattice with free boundary 
conditions. First, we have 

F{tri[L, = 2, = m, free], g) = (g - l)(g - 2)2(— . (15.1) 

For the cases L^^ > 3, we give the result in term of generating functions as shown from eq. 
(9.3) to eq. (9.6). For F{tri[Ly — S^L^ — m, free], q) we calculate 

Af{trt[Ly = 3, free], q, z) = {q - l){q - 2fz[{q - 2) + {q - ifz - {q - ifz^] (15.2) 

V{tri[Ly = 3, free], q, z) ^ 1 - {q^ - 7q^ + 2lq^ - 33g + 23)^ + 2{q - 2fz^ . (15.3) 
There are thus ^ F,tri,Ly='ijree,\ — 2 terms, 

AtH,L,=3,/.eeJ = ^ " V + 21?' - 33? + 23 ± /R^^ ] , J = 1, 2 (15.4) 

where 

R^^^ = g8 _ ^ g^^6 _ 3gQ^5 ^ 949^4 _ yj^^^ ^ 2047?^ - 1486g + 497 . (15.5) 

The locus B for the ^ oo limit of this strip is shown in Fig. 13. This locus is 
compact, consisting of two pairs of complex-conjugate arcs that do not separate the q plane 
into different regions and do not cross the real axis, so that no q^f is defined. The eight 
endpoints of the arcs occur at the eight zeros of -Rts/, where the square root in (15.4) have 
branch point singularities. Thus, we find that the loci B for all of the lattice strips that 
we have considered with free longitudinal boundary conditions (which include the strips 
with free and cylindrical boundary conditions) are compact. These thus contrast with the 
behavior of B for many of the strips with periodic or twisted periodic longitudinal boundary 
conditions, for which B is noncompact. 
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Re(q) 



Figure 13: Singular locus B in the q plane for fl{tri, 3 x oo, free, q) for the 3 x oo strip of the triangular 
lattice with free boundary conditions. For comparison, zeros of the flow polynomial F{tri,3 x L^, free,q) 
for Lx = 21 are also shown. 
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For F{tri[Ly = 4, = m, free], q) we find NF,tri,Ly='ijree,x = 6 and the following coeffi- 
cients of the generating function: 

Atri,4,free,0 = (15.6) 
AtriAfree,! = {q - l)iq - 2f (15.7) 

Atri,4Jree,2 = {q ' l){q ' 2f{2q'^ - 18^^ + 69^^ - UOq^ + 152q^ - 80q + 17) (15.8) 

Atri,4,free,3 = (? - 1) (? - 2)^ (g^O - 1 V + 8 ^ - 33 V + 81 ^ - 1328?^ 

+1520?^ - 131 V + 962g2 - 557^ + 160) (15.9) 

AtriAfreeA = -(? - l)'(? - 2f(q^' - 15?^° + 102?^ - 404g« + 994g^ - 1452?^ 

+856g^ + 990g^ - 2580g2 + 2406^^ - llOOg + 194) (15.10) 

AtriA,free,5 = (? " l)'(? " 2f{2q'' - 37?^° + 319?^ - 1678g« + 5936g^ 
-14675g^ + 25474^5 - 30225?^ + 22685g^ - 8607?^ + 9q + 793) (15.11) 

Atri,4jree,6 = {q ' l)'(g ' 2f{q^^ - 2Aq'^ + 265^1° - 1785g9 + 8181gS - 26897g^ 
+65062g^ - 116663g^ + 153869?^ - 145558g^ + 93785?^ - 37025? 
+6809) (15.12) 

AtriAfreej = (? - lf{q - 2)^(2g' - 19g^ + Ihq^ - 151g' + 155g - 65) (15.13) 

htriAfree,! = -Q^ + 10?^ - 46g^ + 126^^ - 221^2 + 240? - 129 (15.14) 

btri,4jree,2 = - 18?'' + 150?^ - 760?^ + 2592?^ - 6230?^ + 10746g^ 

-13278g3 + 11457g2 - 6404g + 1820 (15.15) 

btriAfree,3 = + 22g^^ - 227?^° + 1448?'' - 6347?^ + 20118g^ - 47275?^ 

+83064g5 - 108597g^ + 103580g2 - 69056g=^ + 29332g - 6141 (15.16) 

btri,4Jree,4 = 2q^^ - Uq^^ + 450g^° - 2824g9 + 12099g^ - 37264g^ + 84596^^ 
-142676?^ + 177551?^ - 159174?=^ + 97739g^ - 36996? + 6557 (15.17) 
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btri,4Jree,5 = - {q - l)\q - 2^ {q^ - IGq' + lUq^ - mq'> + 1225q^ - 2066q^ 

+2201?^ - 1356? + 375) (15.18) 

W/ree,6 = (g - - 2)^ . (15.19) 

We have also calculated the generating function for F{tri[Ly — 5,Lx — m, free], q). The 
results are too lengthy to present here, but we mention that 

deg V{tri, Ly = 5, free) = A^F,tri,5,/ree = 13 . (15.20) 

The respective loci B for the L^, — > oo limits of these strips of the triangular lattice with 
Ly = 4,5 and free boundary conditions again consist of arcs. 

16 Strips of the Triangular Lattice with Cyhndrical 
Boundary Conditions 

We give the result in term of generating functions as shown from eq. (9.3) to eq. (9.6). For 



F{tri[Ly — 2,Lx — m, cyl], q) we have 

Af{trt[Ly = 2, cyl.],q, z) = {q - 1)[1 - {q - l){q - 2){q - ?,)z\ (16.1) 

V{tri[Ly = 2, cyl.], q, z) ^ 1 - {q^ - 6q^ + 15?^ - 19q + n)z - {q - lf{q - 2)z^ . (16.2) 
The coefficients of the generating function for F{tri, Ly — 3,Lx — m, cyl., q) are 

Ah,3,cj//,o = g - 1 (16.3) 

Atri,^,cyi,i = -(? - 1)'(2?" - W + 58g2 - 97g + 68) (16.4) 

Atri,w = -(g - l)\q - 2)(g' - 3g' - 4? + 13) (16.5) 

htrifl,cyi,i = -q^ + - 36g^ + SAq^ - mq^ + 125g - 65 (16.6) 

Wcj/;,2 = {q- l)(g' - llg' + 49g5 - llOg" + 119?^ - 2Aq^ - 69g + 52) (16.7) 

htri,w = - 1)'(? - 2)(g' - 2g - 1) . (16.8) 

We have also calculated the generating function for F{tri[Ly — 4, L^ — m, cyl.], q). The 
results are too lengthy to present here, but we mention that 

deg V(tri, Ly = 4, cyl.) = NF,tri,4,cyi = 6 . (16.9) 
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17 B for 2D Lattices 



We can also obtain some results on Bfi, in particular, qcf values, for (infinite) 2D lattices. 
These follow directly from eqs. (3.4.15). Using the relation (3.4.16) together with the result 
that qdsq) = 3 [80], we have 

Qcfisq) = 3 . (17.1) 
Using (3.4.17) together with the result that qdtri) = 4 [29], we have 

qcfihc) = 4 . (17.2) 

For the honeycomb lattice, formally, qc — (1/2) (3 + y/b) [31, 81], so that 

Q _|_ /E^ 

Icfitri) = = 2.6180... (17.3) 

Since qdkag) — 3, we finally have 

qcfidiced) = 3 . (17.4) 

As was discussed in our previous works on chromatic polynomials, from a study of the 
loci Bw for infinite-length strips of regular lattices of increasing widths (and with a variety 
of boundary conditions) , one can formulate plausible inferences about the the boundary Bw 
for the limit of infinite width, i.e., the full 2D lattice. Given the duality relation (3.4.15), one 
can use either Bw oi Bfi, or both, for this purpose. This program is useful since one does not 
know Bw for any 2D lattice except for results on the triangular lattice ([29], see, however, 
[62]). It was found in [32] and subsequent papers that if one uses periodic longitudinal 
boundary conditions, then Bw exhibits a number of features that one would infer to hold 
for the 2D lattice, such as the property of separating the q plane into various regions, and 
crossing the real q axis at g = and a maximal value, qc- The compactness of the loci 
Bw found for these strips is also the same as the compactness property of Bw for infinite 
regular lattices. This latter property is established by taking the |\/| ^ cxo limit of a section 
of regular lattice A and applying the bound (3.5.1). This compactness is in accord with 
the fact, as discussed in [37, 38], that if Bw did extend infinitely far away from the origin, 
passing through 1/g = 0, then this point would be a point of nonanalyticity of the function 
W{A. q)/q, thereby precluding a Taylor-series expansion of this function around this point. 
However, well-known procedures exist for calculating Taylor-series expansions of W{A,q)/q 
around 1/g = for regular lattices [25, 26]. 

In contrast, as discussed above, we find that the loci Bfi for the infinite-length limits of 
strips with periodic longitudinal boundary conditions are usually noncompact (an exception 
being the case of self-dual strips of the square lattice) and do not, in general, pass through 
q = 0, although they do separate the q plane into various regions. From the duality relation 
(3.4.15) in conjunction with the bound (3.5.1), we infer that Bfi is compact for any infinite 
regular planar lattice A that has a dual planar lattice A* (given the fact that a regular lattice 
has a fixed, finite degree for all of its vertices). In view of this result, there are thus two 
possibilities: (i) either Bfi will become compact for sufficiently great finite width for each 
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type of infinite-length lattice strip, or (ii) wc encounter another kind of noncommutativity in 
addition to (3.4.8), namely that lim^^^oo ^fi{Gs, Ly x oo) is different from the accumulation 
set of the zeros of the flow polynomial of an Ly x section of a regular lattice obtained 
by letting and Ly both approach infinity with Ly/L^ a finite nonzero constant. Further 
study is needed to decide which of these two types of behavior occurs. 

Our present results show that the use of Bw and Bfi calculated on infinite- length finite- 
width strips are somewhat complementary. For example, the study of the loci Bw for finite- 
width, infinite-length strips of a lattice A with periodic longitudinal boundary conditions has 
the advantage that these loci are compact, as is true of Bw{^) = Bfi{\*) for the infinite dual 
pair of lattices A and A*. On the other hand, the loci Bfi for all of the cyclic/Mobius strips 
of the square lattice for which we have calculated them, do exhibit the interesting feature 
of having q^f equal to the value 3 for the infinite square lattice, whereas the Qc values for 
Bw on these strips only approach the square-lattice value asymptotically as Ly increases. Of 
course, for the self-dual families of planar strip graphs studied in [57, 68], Bfi — Bw for each 
value of Ly separately, and in these cases, these loci share the common property of being 
compact and having Qc = Qcf = 3, the asymptotic value. 

18 Summary 

In this paper we have given exact calculations of fiow polynomials F{G, q) for lattice strips of 
various fixed widths and arbitrarily great lengths, with several different boundary conditions. 
We have determined the resultant functions // giving the g-flows per face in the limit of 
infinite-length strips. We have also studied the zeros of F{G, q) in the complex q plane and 
determined exactly the asymptotic accumulation set of these zeros Bfi, across which // is 
nonanalytic in the infinite- length limit. We found that these loci were noncompact for many 
strip graphs with periodic (or twisted periodic) longitudinal boundary conditions, in contrast 
to the usual behavior for the analogous loci Bw for the W function obtained from chromatic 
polynomials for these strips. We also found the interesting feature that, aside from the trivial 
case Ly = 1, the maximal point, qcf, where B crosses the real axis, is universal on cyclic and 
Mobius strips of the square lattice for all widths for which we have calculated it and is equal 
to the asymptotic value qcf — 3 for the infinite square lattice. Duality relations were used 
to derive a number of connections between fl and W, and Bfi and Bw, for planar families 
of graphs. 
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19 Appendix 

19.1 Ly = 2 Cyclic and Mobius Strips of the Square Lattice 

In this appendix we give some explicit flow polynomials for the families of graphs considered 
in the text, in order to illustrate the general discussion in the text. We begin with strips of 



the square lattice. For Ly — 2 we have 

F{sq[2xl,cyc.],q) = Q (19.1.1) 
(this vanishes since this graph contains a bridge) 

F{sq[2 X 2, cyc.],q) = {q - l){q - 2f (19.1.2) 

F{sq[2 X 3, eye], g) = (g - l)(g - 2)(g - 3)^ (19.1.3) 

F{sq[2 X 4, cyclq) = (g - l){q - 2){q^ - 9?" + 29g - 32) (19.1.4) 

F(sq[2 X 5, eye], g) = (g - l)(g - 2)(g - 3)(q^ - 9q^ + 30q - 35) (19.1.5) 

F{sq[2 X 1, Mb.], g) = (g - l)(g - 2) (19.1.6) 

F{sq[2 X 2, Mb.],q) = {q - l)(g - 2)(g - 3) (19.1.7) 

F{sq[2 X 3, Mh.],q) = (g - l)(g - 2)(g^ - 6g + 10) (19.1.8) 

F{sq\2 X 4, Mh.],q) = (g - l)(g - 2)(g - 3)(g2 - 6g + 11) (19.1.9) 

F(sg[2 X 5, Mh.], q) ^ {q - l)(g - 2)(g^ - 12q^ + 57g2 - 125g + 106) (19.1.10) 

19.2 Ly = S Cyclic and Mobius Strips of the Square Lattice 

F(sg[3 X l,q/c.],g) = (19.2.1) 
(where again, this vanishes because the graph contains a bridge) 

F(sg[3 X 2, cyc.],q) = (g - l)(g - 2)\q^ - 3g + 3) (19.2.2) 

F{sq[3 X 3, eye], g) = (g - l)(g - 2)(g - 3f{q' - Qq^ + 14g - 13) (19.2.3) 

F(sg[3 X 4,cyc.],g) = (g - l)(g - 2)(g^ - 17g^ + 129g^ - 569g^ + 1588g=^ 

-2824g2 + 2969g - 1416) (19.2.4) 

F(sg[3 X 5, eye] , g) = (g - 1) (g - 2) (g - 3) (g^ - 19g^ + 165g^ - 860g^ + 2964g^ 

-6977g^ + 11034g2 - 10740g + 4895) (19.2.5) 
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F{sq[3 X 1, Mb.], q) = {q- lf{q - 2) (19.2.6) 
F(sg[3 X 2, M6.], g) = (g - l)(g - 2)(g-^ - Tg^ + 18g - 17) (19.2.7) 
F(sg[3 X 3, M6.], g) = (g - l)(g - 2)(g^ - 12g^ + Glg^ - 167g2 + 253g - 172) (19.2.8) 

F(sg[3 X 4,M6.],g) = (g - l)(g - 2)(g^ - 17g^ + 129g^ - 571g^ 

+1610g^ - 2930g^ + 3222g - 1655) (19.2.9) 

19.3 Ly = ^ Cyclic and Mobius Strips of the Square Lattice 

F(sg[4 X 1, eye], g) = (19.3.1) 

F(sg[4 X 2, cycl q) = {q - l)(g - 2f{q' - 3g + 3)^ (19.3.2) 
F(sg[4 X 3, cyc.],q) = (g - l)(g - 2)(g - Sfiq^ - 6q^ + 14g - 13)^ (19.3.3) 

F(sg[4 X 4, eye], g) = (g - l)(g - 2)(g^^ - 25g^° + 293g^ - 2128g^ + 10664g^ - 38808g^ 
+104899g^ - 211073g^ + 310424g^ - 318191g^ + 204597g - 62416) (19.3.4) 

F(sg[4 X 1, Mb.], g) = (g - l)(g - 2){q'^ - 3g + 3) (19.3.5) 
F(sg[4 X 2, Mb.], g) = (g - l)(g - 2){q^ - llg^ + 52g3 - 133g2 + 188g - 119) (19.3.6) 

F(sg[4 X 3,M6.],g) = (g - l)(g - 2)(g^ - 18g^ + 148g^ - 730g^ + 2380g^ - 5301g^ 

+7968g2 - 7475g + 3369) (19.3.7) 

F(sg[4 X 4, Mb.], g) = (g - l)(g - 2)(g^^ - 25g^° + 293g'' - 2131g^ + 10724g^ - 39374g^ 
+108174g^ - 223727g^ + 343714g3 - 376152g2 + 265280g - 91429) (19.3.8) 

19.4 Ly = 2 Strips of the Square Lattice with Torus and Klein 
Bottle Boundary Conditions 

F(sg[2 X 1, torus], q) = {q - if (19.4.1) 
F(sg[2 X 2, torus], g) = (g - l)(g^ - 7g^ + 21g2 - 29g + 15) (19.4.2) 
F(sg[2 X 3, torus] , g) = (g - 1) (g^ - llg^ + 55g^ - 156g^ + 262g^ - 244g + 97) (19.4.3) 

F(sg[2 X 4, torus] ,q) = (g - 1) (g^ - 15g^ + 105g^ - 441g5 + 1205g^ 
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-2182?^ + 2550g2 - 1755q + 543) (19.4.4) 

F{sq[2 X 1, kb], q) = {q- l)(g^ - 3g + 3) (19.4.5) 
F{sq[2 X 2, kh],q) = {q - l)(g^ - 7q^ + 21^^ - 30g + 17) (19.4.6) 
F{sq[2 X 3, kb], q) ^ {q - l){q^ - llq^ + 55g^ - 156?^ + 262?^ - 245g + 99) (19.4.7) 

F{sq[2x 4:,kb],q) = {q - l){q^ - Aq + 5){q^ - llq^ + 56q^ - 162q^ 

+277?^ - 264g + 109) (19.4.8) 

19.5 Ly = 3 Strips of the Square Lattice with Torus and Klein 
Bottle Boundary Conditions 

F{sq[3 X 1, torus], q) = {q - l)^ (19.5.1) 
F{sq[3 X 2, torus], q) = (q - l)(g^ - llq^ + 55g^ - 156?^ + 262?^ - 244g + 97) (19.5.2) 

F{sq[3 X 3, torus], q) = {q - l){q^ - 17q^ + 136?^ - 671g^ + 2254g^ - 5363g^ + 9076?^ 

-10564^2 + 7655g - 2605) (19.5.3) 

F{sq[3 X 4, torus] ,q) = {q-l) {q^^ - 23q^^ + 253g^° - 1759g^ + 8615g^ - 31363g^ 

+87233?^ - 187346?^ + 309415?^ - 384602?=^ + 342441?^ 
-196170^ + 54481) (19.5.4) 

F{sq[3 X 1, kb],q) = (? - lf{q^ - 3g + 3) (19.5.5) 
F{sq[3 X 2, kb], g) = (? - l){q^ - Hq^ + 55?^ - 158?^ + 275?^ - 273g + 119) (19.5.6) 

F{sq[3 X 3, kb], q) = {q - l){q^ - 17 q^ + 136?^ - 671?^ + 2254g^ - 5365?^ + 9096?^ 

-10645g2 + 7808g-2715) (19.5.7) 

F{sq[3 X 4, kb] ,q) = {q - 1) (g^^ - 23g" + 253g^° - 1759^^ + 8615g® - 31363g^ 

+87233g^ - 187348?^ + 309441g^ - 384752?=^ + 342906?^ 
-196915g + 54959) (19.5.8) 

59 



19.6 Ly = 4 Strips of the Square Lattice with Cyhndrical Boundary 
Conditions 

F{sq[4:Xl,cyl],q)^q-l (19.6.1) 

F{sq[A X 2, cyl],q) = F{sq, 2 x 4, eye, q) ^ (q - l){q - 2)(q' - 9q^ + 29q - 32) (19.6.2) 

F(sg[4 X 3, eyl],q) ^ {q - l){q - 2){q^ - 17 q^ + 129g^ - 569g^ 

+1588?^ - 2824^2 + 2969g - 1416) (19.6.3) 

F{sq[A X 4, cyl\,q) = {q-l){q- 2){q^^ - 25g^° + 293g^ - 2128g* 
+10664?^ - 38808g^ + 104899?^ - 211073g^ + 310424?^ - 318191?^ 

+204597g - 62416) (19.6.4) 

19.7 Ly = 3 Cychc and Mobius Strips of the Honeycomb Lattice 

F{he[3 X 2, cye.],q) = (g - l)(g - 3)(g - 2)=^ (19.7.1) 

F(he[3 X 3, cye.],q) ^ (q - l)(q - 2)(g - 3f{q^ - 9q^ + 29q - 32) (19.7.2) 

F{he[3 X 4, eye.] ,q) = (g - 1) (g - 2) {q - 3) {q^ - 18q^ + 141g^ - 620q^ 

+ 1619g^- 2368g+ 1496) (19.7.3) 

F{he[3 X 5, eye.], q) = {q - l){q -2){q- 3){q^ - 2Aq^ + 260g^ - 1670g^ + 6999g^ 

-19698g^ + 36408g^ - 40240g + 20170) (19.7.4) 

19.8 Ly = 2 Cyclic Strips of the Triangular Lattice 

F{tri[2 X 1, cyc.],q) = (g - 1)^ (19.8.1) 

F{tri[2 X 2, cyc.],q) = (q - l)(g^ - 7g^ + 21q^ - 30g + 17) (19.8.2) 

F{tri[2 X 3, eye.], g) = (g - l)(g^ - llg^ + 55g^ - 159g^ + 282g2 - 290g + 133) (19.8.3) 

F{tri[2 X 4, eye.], g) = (g - l)(g^ - 15g^ + 105g^ - 447g^ + 1269g^ 

-2463g=^ + 3185g^ - 2494g + 897) (19.8.4) 
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19.9 Ly = 3 Cyclic Strips of the Triangular Lattice 

F{tri[3 X 1, cyc.],q) = (q - if (19.9.1) 
F{tri[3x2,cyc.],q) = (g-l)(g^-3g+3)(g^-10g^+45g^-115g^+177g^-157g+63) (19.9.2) 



F{tri[3 X 3, eye], q) = {q - l){q^^ - 2Qq^^ + 190g^° - 113V 

+4743?^ - 14683g^ + 34624?^ - 62969?^ + 88184g^ - 93366?^ + 71377?^ 

-35497? + 8671) (19.9.3) 
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